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The time evolution and the vertical structure of the ocean current is examined 
using moored current meter observations from hurricane Frederic in 1979 and data 
acquired by Airborne expendable Current Profilers (AXCP) in hurricane Norbert in 
1984. Energetic near-inertial oscillations excited by the passage of these hurricanes 
have frequencies that are shifted by 1-20 % above the local inertial frequency. These 
oscillations are evident in the upper 1000 m of the water column and are primarily due 
to the excitation of low-order vertical modes. 

The first five free vertical modr- are calculated from the Brunt-Vaisala frequency 
and the Sturm- Liouville problem. The horizontal velocity eigenfunctions for a flat 
bottom and sloping bottom are fit to the demodulated amplitudes observed in Frederic. 
In the wake of the hurricane, the time evolution of the depth-averaged component plus 
the first two baroclinic modes explain about 60% and 68% of the near-inertial current 
variance using the flat- and sloping-bottom models, respectively. Since the AXCP 
observations were acquired in the high wind regime, the near-inertial response for the 
3-dimensional velocities is simulated by projecting a hurricane-like wind stress field 
onto the first five baroclinic modes. The divergence and curl of the wind stress are also 
convolved with the Green's function (J Q ) for each baroclinic mode. In hurricane 
Norbert, the sum of the first four near-inertial modes explains over 70 % of the 
observed current variance on the right side of the storm track. Most of this current 
variability can be accounted for by the curl terms, although the divergent and Ekman 
terms do contribute to the observed current variance within 30-60 km of the storm 
center. The effect of the stress divergence and Ekman terms on the ocean current 
response rapidly diminishes after one inertial period. 

The role of the depth-averaged velocity component is assessed using numerical 
simulations from a 17-level, primitive equation model with a free-surface boundary 
condition. The slope of the sea-surface induces a depth-averaged velocity that is 
comparable to that observed in Frederic. The simulated modal amplitudes for the first 
two baroclinic modes are within 10-15 % of those derived from the Frederic data. 
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S large scale limit (L/otj'*) 

S t non-dimensional storm translational speed (Price) 

<7 angular wave frequency 
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I. INTRODUCTION 



Thought breeds thought Thoreau, Walden 

As a hurricane moves over the ocean, the wind stress curl causes a divergence of 
upper ocean currents from the storm track, which induces the upwelling of cooler water 
from the thermocline. The magnitude of the ocean's thermal and momentum response 
depends on the scales of the applied atmospheric forcing and the initial ocean 
stratification. Farther back in the wake (rear) of the hurricane, the upper ocean 
currents converge toward the storm track and surface water is downwelled into the 
thermocline. This pattern of upwrilr-g and downwelling regimes in the wake is thus 
part of a time-dependent, three-dimensional response problem. According to linear 
theory, the horizontal wavelength between these cells is proportional to the product of 
the translation speed of the hurricane and the local inertial period. 

The strong atmospheric forcing and the displacement of the isopycnals during the 
upwelling and downwelling phase tend to excite large-scale internal waves in the upper 
ocean. As these waves propagate away from the storm track, the ocean's mass and 
current fields adjust geostrophically to form a ridge along the periphery of the 
upwelling and downwelling regimes. This ridge in the isopycnals and the internal 
waves are part of the ocean's baroclinic response to the passage of a hurricane. The 
wind stress curl also causes a mean mass divergence that depresses the sea-surface 
height relative to the undisturbed height. This depression in the free surface height, 
which is called a barotropic trough, can excite a depth-independent component. 

The forced internal waves are characterized as having large horizontal 
wavelengths, fairly signficant amplitudes, periods close to the inertial period and are 
clockwise rotating. Because the period of the forced waves is so close to the inertial 
period, these oscillations are also called near-inertial. Early theoretical studies 
suggested that the near-inertial response to strong atmospheric forcing should be 
confined to the upper ocean layers (200-300 m). However, recent observations indicate 
that the excitation of the large scale near-inertial waves are associated with low-order 
modes that have vertical wavelengths of 500-1000 m. 

This study focuses on the time-dependent, three-dimensional current response to 
hurricane passage. The emphasis here is on isolating the near-inertial response from 
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the ocean currents and expanding these data into free or forced modes to better 
understand the vertical structure of the ocean response to strong atmospheric forcing. 

A. BACKGROUND 
1. Observations 

One of the first experiments to examine the ocean's thermal response to 
hurricane forcing was during the passage of hurricane Hilda in 1964 (Leipper, 1967). 
Several hydrographic cruises were made to measure the ocean's temperature changes 
using mechanical bathythermographs. The upwelling of cooler water from the 
thermocline was the dominant mechanism along the track, whereas mixing and 
downwelling appeared to have been important outside R max - These measurements 
prompted numerical and analytical investigations of the oceanic response to hurricane 
passage (O'Brien and Reid, 1967; Geisler, 1970). More recently, Black (1983) 
completed an exhaustive 10-year study of ocean temperature observations acquired 
during the passage of hurricanes. Each hurricane was classified according to the storm 
translation speed (U h ), storm intensity (t ), radius of maximum winds (R max ), initial 
ocean stratification and mixed layer depth. One key feature that is relevant to this 
study is that the horizontal wavelength in the thermocline wake was typically 100 - 200 
km instead of 300 - 500 km as predicted by linear theory (Geisler, 1970). Furthermore, 
Black showed that the thermal response directly under the hurricane, which he could 
observe with Airborne expendable BathyThermographs (AXBT), was different from 
that previously observed in the wake. Specifically he found a cresent-shaped pattern of 
maximum temperature decreases of 1-4 °C that extended from the right side of the 
storm to the wake region. There have also been several instances in which the ocean's 
thermal response to hurricane passage has been observed by NOAA data buoys 
(Withee and Johnson, 1976; Johnson and Renwick, 1981). These buoys have measured 
meteorological parameters as well as ocean temperatures, which has contributed to the 
understanding of the upper ocean thermal response. 

Observations of the ocean's current response to the passage of hurricanes have 
been from a few "fortuitous" encounters with current meter moorings that have been 
deployed in support of other experiments. The ocean current speeds have maximum 
amplitudes approaching 100 cm/s and oscillate with frequencies close to the inertial 
period (IP). For example, moored current meter measurements on the continental 
shelf of the mid-Atlantic Bight (100 m) revealed fairly strong current oscillations of 50 - 
70 cm/s (Mayer et al . , 1981). A deeper ocean response was observed during the 
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passage of hurricanes Allen (Brooks, 1983) and Frederic (Shay and Elsberry, 1987a). 
Maximum current amplitudes at 700 m were about 15 cm/s in the Allen measurements, 
whereas the ocean current amplitudes at 950 m in Frederic were 10 cm/s. These sets of 
measurements were acquired on the right side of the storm track where the maximum 
ocean current response is predicted to be the largest (Chang and Anthes, 1978; Price, 
1983). None of the current meter moorings had sufficient vertical resolution to observe 
higher order baroclinic modes in the ocean response to hurricanes. 

Vertical profiles of the ocean current were acquired during the passage of 
midlatitude atmospheric fronts in the Storm Response and Transfer Experiment 
(STREX) using expendable Current Profilers (XCP) (D'Asaro, 1985). Maximum 
ocean currents were about 25 cm/s in the mixed layer, and the current vector rotated 
clockwise with depth in the upper 800 m. This clockwise polarization with depth is 
associated with near-inertial processes and internal waves with large vertical 
wavelengths 0( 1000 m). These large vertical wavelengths are part of the low-mode 
dominance as energy propagates into the interior of the fluid from the wind-forced 
mixed layer. Lai and Sanford (1986) provided a more complete picture of the vertical 
structure of the currents in the far-field after the passage of hurricanes Carrie and 
Dawn in 1972. Moored current meter measurements were augmented with velocity 
profiles in the "Site D" region. They found that over 50 % of the current variability 
could be described using the first three baroclinic modes. 

For the first time, Airborne expendable Current Profilers (AXCP's) were 
deployed within the directly forced or near-field (the positive wind stress curl region) 
during the passage of hurricanes Norbert and Josephine in 1984 and during hurricane 
Gloria in 1985 (Sanford et al, 1987). Not since the early experiments of Leipper (1967) 
has there been such an organized scientific effort to study the ocean response to 
hurricane passage. For the first time, these AXCP's provide a synoptic snapshot of the 
ocean currents with sufficient vertical resolution to adequately describe the hurricane- 
forced baroclinic structure. As will be demonstrated below (see also Shay et al, 1987b) 
over 70 % of the observed current variance in the profiles acquired during hurricane 
Norbert can be explained by a sum of only the first four, forced baroclinic modes. 

2. Analytical and Numerical Modelling 

As the storm translation speed increases, linear theory shows that the ocean 
response is primarily baroclinic and is associated with the dispersion of near-inertial 
waves in the wake (Geisler, 1970). Within the framework of linear dynamics, one of 
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the crucial nondimensional numbers in the oceanic response problem is the Froude 
number, the ratio of the translation speed (U^) to the first baroclinic mode internal 
wave phase speed (cj) (Geisler, 1970). If this Froude number is greater (less) than 
unity, the oceanic response is predominately baroclinic (barotropic). 

In his linear, two-layer treatment of the forced response, Geisler (1970) 
simulated the horizontal structure of the vertical velocity field in the wake of the 
hurricane by convolving a Green's function (Bessel function) with the wind stress curl. 
Upwelling and downwelling regions in the vertical velocity have a wavelength (A) that 
is proportional to the product of the storm translation speed (U h ) and the inertial 
period (IP). Geisler did not explicitly address the structure of the horizontal current 
regime. Geisler points out that it is unlikely that linear theory can describe the ocean 
response in the directly-forced or near-field region. Therefore, extensions of existing 
theory are necessary to facilitate direct comparisons with the AXCP observations 
beneath the hurricane, and to compare the linear theory to simulations from 
sophisticated primitive equation models. 

The wind stress curl induces a semi-permanent baroclinic ridge along the sides 
of the upwelling and downwelling regions. Because this baroclinic ridge may remain in 
the wake for long times following hurricane passage, previous observational and 
theoretical studies have only considered the wind stress curl to be important. However, 
the typical inflow angles in a hurricane are 15-30 °, which induce divergence of the 
wind stress (Holland, 1987). During the passage of hurricane Belle, Mayer et al. (1981) 
demonstrated that the spatial scales of the wind stress curl and divergence were similar 
and that the magnitude of the stress divergence vras about 30 % of that of the wind 
stress curl magnitude. This suggests that the divergence effects in the forcing should be 
examined to see if they contribute to the ocean response beneath the hurricane. 

One of the consequences of linear theory is that the oceanic response can be 
characterized by the first few baroclinic modes in the far- field (Gill, 1984). As 
indicated above, this theoretical prediction appears to be verified by the observations of 
Lai and Sanford (1986) in the far-field of hurricanes Carrie and Dawn. They suggest 
that 75% of the baroclinic, near-inertial energy is contained in the first three bottom- 
slope modes as compared to 55% of the energy in the first three flat-bottom modes. 
Unfortunately, their observations do not measure the barotropic response to hurricane 
passage. The Shay and Elsberry (1987a) study of the near-inertial oscillations observed 
during the passage of hurricane Frederic was limited to a depth-averaged and first two 
baroclinic modes because of the number of current meters in the vertical. 
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Since the baroclinic modes have different rates of horizontal phase velocities, 
it is of interest to study their time evolution. Gill (1984) demonstrated that if the 
storm scale of the wind stress curl exceeds the deformation radius (which is called the 
large-scale limit) of the first baroclinic mode, there will be a time scale over which each 
mode will develop a 71/2 phase difference with the remainder of the solution in the 
mixed layer. This large-scale limit uses a criterion that is not fundamentally different 
from the Froude number of Geisler (1970). A second-order effect is that slight 
increases in the mixed layer kinetic energy will occur on the time scale required for the 
first baroclinic mode to become in-phase with higher order modes. Kundu and 
Thomson (1985) noted this beating effect in numerical simulations of the near-inertial 
response to a moving atmospheric front. As will be shown below (see also Shay and 
Elsberry, 1987a), the time evolution of the near- inertial modal amplitudes from 
hurricane Frederic agrees rather well with these baroclinic time scales for modes 1 and 
2. Moreover, amplitudes of the first two baroclinic modes in the mixed layer are within 
10-15 % of those expected from linear theory. 

The nonlinear, baroclinic response in the ocean to a moving hurricane was 
numerically simulated by Chang and Anthes (1978), Adamec et al. (1981), Price (1981, 
1983), Hopkins (1982), and Greatbatch (1983). In these studies, a rigid lid was 
imposed at the sea surface, which eliminated the barotropic or depth-averaged 
response. For example, Hopkins (1982) simulated the dynamic response to a hurricane 
moving at the same translation speed as Frederic and found that the simulated mixed 
layer currents agreed well with the Shay (1983) observations. However, the computed 
magnitudes of the currents in the thermocline were too small. It is believed that one 
source of the discrepancies between the simulations and the observations is the neglect 
of the depth-averaged mode by imposing a rigid lid in the model. Although Brooks 
(1983) noted a significant depth-averaged current in the wake of hurricane Allen, the 
contribution of the depth-averaged component to the total observed variance was not 
estimated. It will also be demonstrated below (following Shay and Elsberry, 1987a) 
that the near-inertial response to hurricane Frederic contains a depth-independent 
component of 7-10 cm/s that contributes to about 15-20 % of the near-inertial 
variability. 

Both the ocean's baroclinic and barotropic response to hurricane forcing was 
numerically simulated for the first time by Chang (1985). An inherent problem of free- 
surface modelling is that the barotropic mode propagates so quickly that very small 
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model time steps must be used. To circumvent this computational problem, Chang 
(1985) introduced a mode-splitting technique into the free-surface, primitive equation 
model. The simulated current response extended from the surface to 2000 m with a 
barotropic component of about 10-15 cm/s. New numerical model experiments with 
the Chang model will be used in this study. 

B. OBJECTIVE AND KEY ISSUES 

The objective of this study is to understand the vertical structure of the near- 
inertial ocean current response to the passage of hurricanes. The emphasis here is to 
blend observations from hurricanes Frederic and Norbert with analytical and numerical 
models. Thus, the near-inertial response is studied in terms of kinematics, energetics, 
normal modes and wavelengths using free modes and forced models in the far-field and 
near-field, respectively. 

A synopsis of the moored ocean current observations acquired in hurricane 
Frederic in 1979 and current profiles from AXCP's deployed in hurricane Norbert is ' 
given in Chapter II. Brief descriptions of the hurricanes are also given in terms of the 
scales of the applied atmospheric forcing. 

In Chapter III, the near-inertial response to a moving hurricane is examined 
within the framework of linear theory by extending the two-layer model of Geisler 
(1970) to a continuously stratified fluid and superposing the wind stress onto the first 
five baroclinic modes as in Kundu and Thomson (1985). In this treatment, analytical 
expressions are derived for the horizontal current velocities as well as the vertical 
velocity for the "near-field". The effects of the wind stress and the wind stress 
divergence are also included in the analysis. 

The time-evolution of the near-inertial response to hurricane Frederic is examined 
using the moored current meter arrays in Chapter IV. The behavior of the vertical 
structure of the first two baroclinic modes is compared to predictions from linear 
theory (Gill, 1984). The effect of bottom topography on the near-inertial current 
response is also explored using the Lai and Sanford (1986) model. 

The complete vertical structure of the near-inertial response in the near-field is 
addressed in Chapter V using the vertical profiles of ocean currents and temperatures 
acquired in hurricane Norbert. The near-inertial current profiles are expanded in terms 
of the forced dynamical modes using the theories derived in Chapter III. The effect of 
the stress divergence and wind stress are shown to be important in the direct forcing 
region as well as the wind stress curl. 



23 



A controversial issue raised by the Frederic observations (Shay and Elsberrv, 
1987a) is the role of the depth-averaged circulations induced by the passage of 
hurricanes. The role of the depth-averaged current and the baroclinic structure is 
assessed in Chapter VI using the numerical simulations from a 17-level, primitive 
equation model (Chang, 1985). Since the model current simulations appear to have 
sufficient veracity, the model simulations are compared to the Frederic observations. 
Finally, concluding statements are given in Chapter VII. 
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II. DATA DESCRIPTION 



A. HURRICANE FREDERIC 

1. Storm Track and Wind Field 

The storm track shown in Fig. 2.1 is based on the post-season analysis of 
hurricane Frederic's movement across the Gulf of Mexico (Hebert, 1979). Because of 
the intensity of Frederic, reconnaissance aircraft constantly monitored the storm. 
Hurricane Frederic reached maximum wind speeds within 80 to 130 km west of the 
array sites about 21 Universal Time Coordinate (UTC) 12 September (the time of 
closest approach). Visible satellite imagery at 20 UTC 12 September clearly delineated 
a well-developed eye of 50 to 60 km in diameter. The translational speed of the 
hurricane (U^) as it approached the Gulf Coast was 6-7 m/s, which was larger than 
the internal wave phase speed of the first baroclinic mode (cj = '3 m/s). 




Figure 2.1 The path (circles) of hurricane Frederic in the northern Gulf of Mexico 
(Hebert, 1979) beginning 12 UTC 11 Sept, and ending 03 UTC 13 Sept. 1979. The 
boxes depict the XOAA data buoy 42003, the OTEC and three NAVOCEANO (CMA 
1,2 and 3) current meter arrays. 
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Powell (1982) documented the boundary layer wind field associated with 
hurricane Frederic using data collected by two NOAA aircraft. The wind speeds 
ranged between 48 and 58 m/s at a radius of maximum wind (R max ) between 27 and 
33 km. The curl of the wind stress field deduced from these values has a horizontal 
scale of about 60 km. Other wind reports confirm these intensities as Frederic made 
landfall (Hebert, 1979). 

Marine winds (Fig. 2.2) were also measured by NOAA Buoy 42003 (Johnson 

and Renwick, 1981). The eye clearly passed over the buoy, as indicated by the 

minimum in the wind field as the direction changed from 40 to 200 °T, with a 

corresponding decrease in pressure to 959 millibars (mb). The maximum (near-surface) 

wind speeds at this location never exceeded 35 m/s. Using a Rankine vortex and a 

maximum wind stress of 3.5 N/m 2 in Frederic (Black, 1983), the maximum curl is 220 

lO*^ cm / s 2 w j t k a 2R max scale of 60 km. For an inflow angle of 20 °, the maximum 

7 7 

stress convergence is about 75 x 10 cm/s . 

2. Current Meter Data and Measurement Errors 

The local variations in the bottom topography of the DeSoto Canyon (Fig. 
2.3) are approximately parallel and normal to the coast at Current Meter Array 
(CMA)2 and CMA3. For example, north-south flow at CMA2 is in the cross-shelf 
direction, but it is in the along-shelf direction at CMA3. Since the storm moved on a 
northward track and the isobaths are nearly east-west and north-south along the 
northern and eastern rim of the DeSoto Canyon, the ocean current data are 
represented in a standard Cartesian coordinate system (i.e., y positive north and x 
positive east). 

Seventeen Aanderaa RCM-5 current meters were deployed on five moored, 
taut- wire, subsurface arrays 1 (Fig. 2.3) in depths ranging from 100 - 500 m in the 
northern Gulf of Mexico. The arrays were designed with a flotation device one meter 
above the near-surface current meter (nominally 20 m). These current meters sampled 
ocean current, direction and temperature at 10 minute intervals. Two current meter 
arrays (CMA2 and CMA3) were deployed on adjacent sides of the DeSoto Canyon. 
The other NAVOCEANO mooring was on the continental shelf in about 100 m of 
water. 



1 Since the positions (but not the observations) of CMA6 and CMA9 remain 
classified confidential, the results from these arrays arc only compared to the results 

from the unclassified arrays CMA1, CMA2 and CMA3. 
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JULIAN DAY/TIME (GMT) 



Figure 2.2 Wind speed (upper) and direction (lower) measured by NOAA buoy 42003 
from 12 UTC 11 Sept, to 15 UTC 13 Sept. 1979. The hatched area depicts the period 
when the data buoy was in the eye of the hurricane (Johnson and Renwick, 1981). 

The OTEC buoy was located 210 km south-southeast of the DeSoto Canyon 
in 1050 m of water and was instrumented with four Aanderaa RCM-5 current meters 
that measured ocean current speed, direction, temperature and pressure at 20 minute 
intervals (Starr and Maul, 1982). Table 1 is a synopsis of the observations for the 
entire period of deployment, which extends from the end of July to mid-October. The 
storm period is from a few days prior to hurricane passage to the end of the 
deployment period (mid-October). 
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Figure 2.3 Positions of NAVOCEANO current meter arrays (triangles) and AXBT 
observations 1,11 and 12 (rectangles) of Black (1983) relative to the track of hurricane 
Frederic. The minimum pressure and maximum wind speeds at hourly intervals are 
taken from Hebert (1979). The radius of the maximum winds is indicated by the dashed 
lines parallel to the storm track and the depth contours (fathoms) are shown as solid 
lines. 



Because the Savonius rotors were eventually lost from all the current meters in 
the mixed layer due to the large current speeds, the lengths of the time series are not 
equal. Furthermore, these large current speeds have not been corrected for rotor 
pumping. Current speeds recorded by Savonius type meters can be altered by 
horizontal and vertical mooring motion induced by the surface gravity wave field. 

Other sources of errors in the Aanderaa current meters are the instantaneous 
directional sampling, lack of vector averaging, and the slow vane response that 
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contaminates the spectrum of wave motions, especially at the higher frequencies 
(Beardsley et al., 1977). 



TABLE 1 

A synopsis of moored current meter observations 
in hurricane Frederic based on the Universal Time Coordinate (UTC). 



Meter 


Record 


Start 


End 




Depth 
(m j 


Length 

(days) 


Time 

(UTC) 


Time 

(UTC) 


Variables 


CMA1 


21 


47.3 


1600 26 Jul. 


0000 12 Sep. 


u,v 


49 


62.5 


1600 26 Jul. 


0600 27 Sep. 


u,v 


64 


57.1 


1600 26 Jul. 


1920 21 Sep. 


u,v,T 


92* 


61.2 


1600 26 Jul. 


2040 25 Sep. 


u,v,T 


CMA2 


19 


64.9 


1800 26 Jul. 


1650 29 Sep. 


u,v 


179 


65.3 


1800 26 Jul. 


0100 30 Sep. 


u,v,T 


324 


66.0 


1800 26 Jul. 


1900 29 Sep. 


u!v,T 


CMA3 


21 


57.3 


1330 31 Jul. 


2200 26 Sep. 


u,v 


251 


68.5 


1330 31 Jul. 


0100 8 Oct. 


u,v,T 


437 


68.2 


1330 31 Jul. 


1730 7 Oct. 


u|v,T 


457 


64.4 


1330 31 Jul. 


0000 4 Oct. 


u,v,T 


CMA6 


20 


64.9 


2000 26 Jul. 


1715 29 Sep. 


u,v 


180 


61.3 


2000 26 Jul. 


0410 26 Sep. 


u,v,T 


330 


65.6 


2000 26 Jul. 


1150 1 Oct. 


u,v,T 


CMA9 


20 


60.3 


1640 31 Jul. 


0000 29 Sep. 


u,v 


250 


71.4 


1640 31 Jul. 


0320 11 Oct. 


u,v,T 


500 


64.2 


1640 31 Jul. 


2250 3 Oct. 


u,v,T 


OTEC 


100 


30.00 


0007 1 Sep. 


0007 1 Oct. 


u,v,T,p 


232 


30.00 


0007 1 Sep. 


0007 1 Oct. 


u.v,T,p 


546 


30.00 


0007 1 Sep. 


0007 1 Oct. 


u,v T,p 


950 30.00 0007 1 Sep. 

* time clock synchronization problems 


0007 1 Oct. 


u,v,T,p 
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The pressure records on the OTEC current meters indicate that the mooring 
was tilted by the hurricane-forced currents. The maximum tilt of the mooring based on 
the vertical displacements measured by the pressure sensors was about 15 °. The 
current meters at 100, 232, 546 and 950 m were displaced downward by 65, 55, 30 and 
0 m, respectively. The corresponding changes in temperature were about 2.6, 1.5, 0.75 
and 0.25 0 C. 

On the NAVOCEANO arrays, the temperature changes due to vertical 
displacements are not known as pressure sensors were not included. An upper bound 
on the mooring motion can be estimated from the AXBT and moored temperature 
measurements. The NAVOCEANO arrays were designed to withstand currents in 
excess of 150 cm/s with a maximum tilt of 10° (C. Robinson, 1986, personal 
communication). The maximum vertical displacement in the 10° tilt case would be 14 
m and cause the maximum depth of the upper current meter to be 35 m. The AXBT 
observations (see Fig. 2.3) indicated that the threshold temperature of the thermistors 
(21.5 0 C) was located between 75 to 80 m depth. Because the mixed layer (initial 
depth of 21 m) temperatures were always above this threshold, the mixed layer current 
meters must have always remained above 80 m. A 60 m displacement of this 
instrument would correspond to a tilt of the moored array by 21.5 0 from the vertical. 
The corresponding maximum displacement is about 160 m in the horizontal. Thus, the 
horizontal displacements of the mooring array accounts for very little of the observed 
temperature variations. 

Since the mixed layer current measurements were made at 2 1 m below the 
surface, the forced surface gravity wave field will induce a current that will contaminate 
the current measurements (Halpem et al., 1974). The magnitude of the surface wave 
current is estimated from 

U s = (nH s /L s ) 2 c 0 e k s z , (2.1) 

where H s is the significant wave height, L s is the dominant wavelength of the surface 
gravity wave field k s is the significant wavenumber (2 tt/L s ) and c Q is the deep water, 
surface gravity wave phase speed (g/k $ )* The wavelength is estimated from the wave 
period (T s ) using the relationship (g/27t)T s . The significant wave height and period 
of the surface gravity waves were estimated at NOAA data buoy 42003 (Johnson and 
Renwick, 1981). The significant wave height at the buoy was approximately 8 m with 
a period of 10 s . The current induced by the surface wave is estimated to be about 40 
cm/s, which then decreases exponentially with an e-folding depth of 25 m. Since the 
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mooring may have been tilting a maximum of 15 °, the magnitude of the surface wave 
current ranged from 18 and 6 cm/s for tilts of the mooring from 0 - 15 °, respectively. 

The time-dependent behavior of the currents induced by surface waves was 
studied theoretically by Hasselman (1970) and Pollard (1970). As the surface waves 
increase over the first quarter of an inertial period, the current begins to rotate 
inertially due to the influence of the Coriolis force. During the passage of storms (time 
scales of 1/0* surface waves can induce fairly large inertial oscillations (Hasselman, 
1970). The estimated surface wave current component was removed from the mixed 
layer record by forming a time series of an inertially rotating current with amplitudes of 
10 cm/s and phase of n/2 starting 5 h before the time of closest approach. The 
adjustment for the surface wave current was e-folded over 4 IP's following the 
maximum value at CMA3. The phase of the forced near-inertial waves was preserved 
after the removal of this current. 

3. Temperature Profiles 

The AXBT data were collected by Black (1983) during and subsequent to the 
passage of Frederic in the area of the DeSoto Canyon (Fig. 2.3). Since the AXBT's 
only extended to 200 m, it was necessary to extrapolate the vertical temperature 
profiles to the bottom using the moored temperature data. For example, the 
temperature observations at CMA2 were used to extrapolate the thermal profile at 
AXBT station 1 (Fig. 2.4a). The vertical temperature gradients below 200 m were also 
checked against CTD data collected in support of OTEC during the summer of 1977 
(Thomas et al., 1979) to insure consistency. 

The mixed layer temperatures were not available because the ocean 
temperatures exceeded the maximum temperature resolvable by the thermistors. The 
AXBT observations of Black (1983) indicate that the mixed layer cooled by 3°C and 
deepened to 40-50 m in the DeSoto Canyon region. The mixed layer temperatures 
after storm passage were 25-26 °C. The vertical profiles of the Brunt-Vaisala frequency 
(N) were estimated using the AXBT data and climatological temperature-salinity (T-S) 
relationships (Wahl and Teague, 1983). A time-averaged value (Fig. 2.4b) was 
computed at each depth using the N profiles. 

B. HURRICANE NORBERT 

1. Storm Track and Wind Field 

During the summer of 1984, the Hurricane Research Division (HRD) of 
NOAA conducted Planetary Boundary Layer Experiments in hurricane Norbert. Prior 
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to the AXCP research flights (Sanford ei al., 19S7), tropical storm Norbert was slowly 
developing into a hurricane (Fig. 2.5). Norbert slowly drifted in a cyclonic path for 
about 6 days. On 21 September, tropical storm Norbert was classified as a hurricane 
with a minimum central pressure of 950 mb and maximum winds of 1 15 knots. Notice 
the recurving path of hurricane Norbert between the 21 and 23 September just before 
the AXCP research flights on 23 and 24 September. The translation speed (U h ) for 
hurricane Norbert was only 3.8 - 4.2 m/s in the vicinity of the AXCP research flight 

and the vortex was moving at 320°T. The AXCP portion of the experiment w^as 

✓ 

conducted by Horizon Marine with the support of a consortium of oil companies under 
the project Ocean Response to Hurricane Joint Industry Program (ORHJIP). The 
locations of the 16 AXCP's that were successfully deployed in hurricane Norbert are 
shown in Fig. 2.6. 




LONGITUDE 

Figure 2.5 Track of hurricane Norbert and the location of the AXCP research flights 
(from Sanford et al. 1987). The Norbert AXCP flight started at 2230 UTC 23 
September and ended at 0130 UTC 24 September 1984. 
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Figure 2.6 Positions of AXCP's (probe number is indicated) relative to storm center 
of Norbert (19°N, 109 °W, moving towards 320°T). The spatial distribution of the 
AXCP's is scaled nondimensionally in terms of the primary radii of maximum winds 
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There are primary and secondary wind maxima that lead to a double-eye type 
structure as shown in Fig. 2.7 (Willoughby et ai, 1984). The winds in Norbert 
exceeded 50 m/s at both the primary (34 km) and secondary (17 km) radii of maximum 
wind, and the central pressure was in the 950 mb range. Outside the secondary 
maximum, the wind decreases to 20 m/s over 120 km. The maximum wind stress 
estimated for Norbert is 4.4 N/m 2 at the primary radius of maximum winds. 

2. Model of the Wind Stress 

Because of the lack of wind profiles in a few of the quadrants and the double- 
eye structure that would add considerable complexity in the wind field, a simpler 
distribution of the wind is treated here. A common representation of a hurricane wind 
stress distribution is the Rankine vortex given by Chang and Anthes (1978). The 
tangential and radial wind stress components are given by 
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Figure 2.7 Flight-level wind profiles from hurricane Norbert (solid) in the right-rear 
(1), right (2) and (3) right-front quadrants. 

where R max (radius of maximum wind) = 30 km and R out (radius to the outer edge 
of the hurricane) = 300 km. Typical values for |t.| and ItqI used in the Chang and 
Anthes experiments were 1.0 and 3.0 Nt/m 2 , respectively. The function f(0) represents 
the asymmetric factor given by 

f{0) = 1.0 + b cos (0) , 

where b = 0 for the symmetric and b = 0.3 for the asymmetric case. This wind stress 
distribution is input as a body force in the mixed layer. The effect of the asymmetric 
storm was to displace the maximum ocean currents and temperature perturbations 
even farther to the right of the track. This rightward bias produced by the asymmetric 
storm was small compared to the bias induced by nonlinearities in the Chang and 
Anthes model. In this research, the storm is assumed to be symmetric (b = 0). The 
region of strong positive vorticity (wind stress curl) is confined within ± R max and 
corresponds to a length scale of 2R max (68 km for Norbert and 60 km for Frederic). 

The wind stress curl exceeds 280 x 10'^ cm's 2 in the core of the hurricane 
(Fig. 2.8). The near-field is defined as the region where the wind stress curl is greater 
than zero which is approximately at a r fl = 3R max . Since both hurricanes had inflow 
angles of about 20 °, the maximum divergence of the wind stress is about 50 x 10 
cm's 2 . Even larger inflow angles may be needed for the growth and maintenance of 
the tropical cyclone as it moves over the ocean (Holland, 1987). If the inflow angle 
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was 30 °, the maximum divergence would be 100 x 10 cm/s 2 . Thus, the divergence 
term within a few radii of maximum winds may be important, as noted by Mayer et al. 
(1981) during the passage of hurricane Belle. The length scales associated with the 
divergence are comparable to those associated with the wind stress curl (2 R max )- 

3. Ocean Current Data 

A total of 16 AXCP's were sucessfully deployed in hurricane Norbert from the 
NOAA aircraft (Table 2). The profilers sample relative currents once every revolution 
(about 0.3 m in the vertical) with rms errors in the velocity measurements of about 1 
cm/s (Sanford et al., 1982). Since the profiler measures vertical current shear, the 
currents are integrated upward from a reference value as described in Sanford et al. 
(1987). The fast probes descend at a rate of 4.5 m/s and sample currents and 
temperatures in the upper 1000-1500 m of the ocean, whereas the slower probes (2.2 
m/s) provide excellent resolution of the ocean response in the upper ocean (typically 
200-250 m). 

Upper ocean current measurements are contaminated by storm-generated 
surface waves that induce a current. Sanford et al. (1987) developed a least-squares 
model to account for the currents associated with surface wave, which for the u- 
component is 

u m =[A u cos(<Tt) + B u sin(<Tt)]e kz , (2.3) 

where A u and B u are the least squares coefficients (Table 3) for the u-velocity 
component, <T is the surface wave frequency (27t/T), T is wave period (nominally 5-15 
s), k is the wavenumber (cr 2 /g) and g is the acceleration of gravity. The least-squares 
model for the v-component of velocity is similar to the above expression. Equation 
(2.3) is analogous to (2.1) except that the amplitudes of the surface waves (H s ,L s ) are 
replaced by the amplitudes of the surface wave currents (A u ,B u ). Although the 
dominant waves appear to fit the model quite well (Fig. 2.9), some surface wave signals 
may remain in the profiles after removal of the surface wave currents, which may 
possibly contaminate the vertical structure analyses below. However, a second fit to 
the profiles after the dominant surface wave is removed reveals very small amplitudes 
in most of the profiles. 

The velocity profiles are low-pass filtered every 3 m by a simple running mean 
(10 values) for subsequent processing. The AXCP profiles contain vertical structure 
that has not been previously available from moored current meter arrays (Mayer et al., 
1981; Brooks, 1983; Shay and Elsberry, 1987a). For the first time, the vertical structure 
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for hurricane Norbert based on a Rankine vortex and observed storm parameters. The 
abscissa and ordinate axes are scaled in terms of R relative to the storm center as 



TABLE 2 



A synopsis of AXCP observations acquired during hurricane 
Norbert for the fast (F:4.5 m/s) and slow (S: 2.2 m/s) falling probes. 



Probe 


Time (JD) 


Lat.(N) 


Long.(W) 


W s (m/s) 


Rate 


N02 


266.947 


20 34 


108 12 


25.3 


F 


N04 


266.950 


20 15 


108 08 


32.0 


F 


N06 


266.966 


19 36 


108 20 


35.2 


F 


N13 


266.983 


19 04 


108 19 


34.4 


F 


N14 


266.985 


18 58 


108 04 


28.1 


F 


N16 


267.000 


19 12 


108 57 


33.0 


S 


N18 


267.004 


18 50 


108 53 


40.2 


F 


N20 


267.007 


18 31 


108 50 


29.7 


F 


N21 


267.009 


18 20 


108 48 


24.0 


F 


N03 


267.012 


18 07 


108 46 


21.8 


F 


N15 


267.021 


19 25 


109 05 


5.7 


S 


N22 


267.021 


18 53 


108 56 


41.9 


S 


N23 


267.022 


18 54 


108 56 


42.7 


S 


N24 


267.022 


18 55 


108 56 


42.9 


S 


N26 


267.041 


18 56 


109 41 


30.2 


F 


N31 


267.058 


19 49 


109 24 


38.0 


F 



of the ocean current response to hurricanes can be examined in detail with no vertical 
aliasing. 

4. Ocean Temperature Data 

The ocean temperatures measured by the AXCP's are low-pass filtered in a 
similar manner to the velocity data (Fig. 2.10). Mixed layer depths range between 18 
to more than 40 m within the domain of hurricane Norbert. The temperature profile 
from AXCP 20 indicates a large temperature decrease of 1 1 °C in the upper 100 m of 
the water column. This strong gradient impedes the vertical penetration of the 
hurricane-induced response by limiting the entrainment velocity. Vertical scales 
associated with the temperature gradients in the upper thermocline are about 160 m, 
below which the temperature decreases uniformly with depth in the domain influenced 
by hurricane Norbert. 
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Figure 2.9 Original (solid), surface wave (dashed) and residual (dotted) v-component 
current (cm/s) profile from 20 to 200 m for AXCP 04 in Norbert. The residual current 
is obtained by removing the surface wave component from the observed current profile, 
which begins at about 21 m depth. 

The vertical profiles of the Brunt- Vaisalii frequency in Fig. 2.11 are estimated 
from the temperature data using climatological temperature-salinity (T-S) relationships, 
which may induce an error in the Brunt- Vaisala frequency of 0.-0.2 cph over a depth 
range of 100-450 m (Wahl and Teague, 1983). Only a few temperature and salinity 
profiles are available to estimate density structure in the eastern North Pacific where 
the AXCP's were deployed. The spatially-averaged Brunt-Vaisala frequency increased 
to about 8 cph at the base of the mixed layer and decreased to 2 cph at 200 m and 
below. The stratification suggests that the vertical penetration of the ocean's thermal 
response induced by the hurricane may be limited to the upper ocean. 
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TABLE 3 



Periods and velocity amplitudes for the u 
(A u ,B u ) and v (A V ,B V ) components of the surface wave in hurricane Norbert. 
based on the model of (Sanford et ai, 19S7). 



Probe Period A u 





(s) 


(cm/s) 


N02 


9 


• 24 


N04 


10 


24- 


N06 


9 


-95 


N13 


9 


-51 


N14 


10 


-70 


N15 


10 


-25 


N16 


11 


-19 


N18 


10 


-129 


N20 


10 


-50 


N21 


10 


-37 


N03 


9 


-25 


N23 


10 


-57 


N24 


10 


-67 


N26 


9 


-147 


N31 


9 


0 



Bu 


\ 


B v 


(cm/s) 


(cm/s) 


(cm/s) 


11 


4 


-61 


-35 


-95 


122 


26 


113 


196 


-39 


-26 


-43 


-50 


12 


-58 


-13 


-5 


5 


-5 


8 


-29 


1 


20 


-51 


43 


50 


81 


-31 


-123 


-13 


-76 


-21 


38 


13 


65 


28 


-1 


9 


-46 


117 


-68 


-60 


107 


5 


84 



C. AIR/SEA PARAMETERS 

The initial oceanic response (spin-up) to hurricanes is governed by the 
parameters (Table 4) of the applied atmospheric forcing (Geisler, 1970; Price, 1984; 
Greatbatch, 1984). The translation speeds (U h ) for hurricanes Norbert and Frederic 
are only 4 - 6.5 m/s. Hurricane Norbert is characterized as a slowly moving storm. 

The maximum wind stress estimated for Norbert is 4.4 N/m 2 at the primary radius of 
maximum winds. By contrast, Frederic is a fast moving storm at 6.5 in's, but it is less 
intense than Norbert with a t of 3.5 N/m 2 . 

1. Ocean Parameters and Nondimensional Framework 

The initial mixed layer depths (h) in the region of hurricanes Norbert and 
Frederic were 30 to 45 m (see Table 4), and thicknesses of the thermocline (b) were 160 
and 200 m, respectively. Both storms were located in the subtropics and had /nertial 
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Figure 2.10 Low-pass filtered temperature (°C) profile from hurricane Norbert AXCP 
20 . 

Period's (IP) ranging from a day (Frederic) to a day and a half (Norbert). The first 
mode internal wave phase speeds are 2.1 and 3 m/s based on the spatially-averaged 
Brunt-Vaisala frequency and the Sturm- Liouville problem (Chap. III). 

The important nondimensional numbers in the ocean spin-up to hurricane 
passage are associated with comparisons between atmospheric and oceanic scales. 

Geisler (1970) demonstrated that if the internal Froude number (C T h /Cj) > ( < ) 1, the 
oceanic response is primarily baroclinic (barotropic). Within this framework, a mixed 
baroclinic and barotropic response may be expected since the Froude number is about 
1.9 in Norbert and 2.2 in Frederic. Another relevant scale intercomparison is the ratio 
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Figure 2.11 Spatially-averaged Brunt-Vaisala frequency (cph) profile from Norbert with 
the 95% confidence intervals at selected depths based on the bootstrap method and a 
normal distribution. 

of the scale of the wind stress curl (2R max ) to the deformation radius of the first 
baroclinic mode (a^ (Cj/f). If this nondimensional forcing scale is > > (<) 1, the 
response is primarily baroclinic (barotropic). As with the internal Froude number, a 
mixed baroclinic-barotropic response is expected in these cases. Although the internal 
Froude number and nondimensional forcing scale are not fundamentally different, the 
nondimensional forcing scale is also the large-scale limit for which linear theory should 
apply (Gill, 1984). Veronis (1956) showed that if the ratio of the nondimensional 
forcing scale to the Froude number is larger than 7t, the internal wave amplitude 
associated with the baroclinic response will be large. Because this ratio 
(nondimensional time scale) is 0(1) in both of these storms, the amplitude of the 
internal waves is smaller in Norbert because Frederic is considered to be a fast moving 
storm. 

The Geisler (1970) wavelength (A) is shown to be a function of storm 
translational speed, the Coriolis parameter, and the phase speed of the first baroclinic 
mode (cj): 
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TABLE 4 

Air/sea parameters for hurricanes Norbert and Frederic 



^max 


Norbert 

34 


Frederic 

30 


T (Nt/m 2 ) 


4.4 


3.5 


U h (m/s) 


4.0 


6.5 


h (m) 


40 ' 


45 


g' x lO'-qm/s) 


3.7 


2.2 


b (m) 


160 


200 


IP (d) 


1.5 


1.0 


Cj (m/s) 


2.1 


3. 


A (km) 


444 


520 


U h/fl <F) 

2 W«r' (?) 


1.9 


2.2 


44 


41 


1.5 


1.5 


(S)/(F) 


0.9 


1.7 



A = WfXV-c, 2 ) 1 / 2 . 

According to the above expression, as Cj approaches the wavelength of the near- 
inertial waves decreases. 

Price (1984) discussed the importance of the nondimensional numbers derived 
from the air-sea parameters. An important nondimensional number in the mixed layer 
is the Burger number M, which is a measure of the importance of the horizontal 
pressure gradients 



M = 



(I + l/S, 2 )g -h 

( 2R max0 2 ’ 



(2.4) 



where S t is the nondimensional storm speed (U h /2R max f), h is the mixed layer depth 
and g' is the reduced gravity. The Burger number measures the importance of the 
horizontal pressure gradient that couples the mixed layer to the thermocline and causes 
a vertical transfer of energy from the wind-forced mixed layer into the thermocline. 



Price (1984) suggests that the shift in the frequency above the local inertial frequency 
(blue-shift) is equal to one half of the mixed layer Burger number or Av. The Burger 
numbers in Frederic and Norbert are approximately 0.08 and 0.36. The associated 
frequency shifts of the near-inertial response should be 0.04 and 0.18 above the local 
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inertial frequency. These shifts are within the near-inertial frequency band of 0.9f-l.2f 
(Kunze, 1985; Mooers, 1975). Thus, the near-inertial frequencies will be set to 1.18f 
for hurricane Norbert in light of the excellent agreement between the estimated Burger 
numbers and the frequency of the near-inertial oscillations in hurricane Frederic (Shay 
and Elsberry, 1987a). 

Price (1984) also defined two other important nondimensional parameters in 
the response problem: 

S = — and 

^max^ 

Mb 

T= r* 

h l 

where b is th-: thermocline thickness, and g' is the reduced gravity (Table 4). These 
non-dimensional numbers are given in Table 5, using the air- sea parameters specified in 
Table 4. The velocity in the thermocline depends on the nondimensional storm speed, 
S t< For a symmetric storm, S t is a ratio of the inertial time scale f * to the advective 
time scale of the storm, 2R max /Uj 1 . Alternatively, S t may be intrepreted as a ratio of 
the along-track wavelength U^/f to the scale of the wind stress curl (2R max ). The 
parameter S t is also the nondimensional number k f derived by Greatbatch (1984) who 
refers to k. as the translational storm scale (ratio of the local inertial period to the time 
scale over which the forcing acts). The translational time scale can also be thought of 
as the time available for mixing. 

A storm with S t < 1 is considered to be slow whereas a fast storm has S t > 2. 
According to the Price criterion, Norbert would be considered a fast storm (S t = 2.1), 
whereas Frederic is a transition storm since S t = 1.7. The value of S t also sets the 
amplitude ratio between the geostrophic and wave components associated with the 
near-inertial response. According to Black's (1983) classification of storms, Norbert 
would be considered a fast moving storm (U^> 3 m/s). If the advective time scale is 
long compared to the inertial time scale (S t < 1, i.e. a slow storm), the atmospheric 
forcing can decrease the amount of kinetic energy in the mixed layer in the rear half of 
the storm (Price, 1984). 

The thermocline Burger number T represents the vertical phase difference of 
the velocity across the thermocline. This vertical phase difference provides the 
mechanism for the vertical propagation of energy from the wind-forced mixed layer to 
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TABLE 5 



Nondimensional numbers from Price (1984) model 
versus hurricanes Frederic and Norbert. 



# 


Price 


Frederic 


Norbert 


S t 


1.42 


1.74 


2.07 


M 


0.12 


0.08 


0.36 


T 


0.48 


0.27 


2.10 


R* 


0.28 


0.18 


0.47 


A 


0.01 


0.1 


0.01 



A is the vertical aspect ratio 
R is the Rossby number 

* These values were computed using the expression T/hjU^f 
because Price (1983) includes a constant of 3 in the numerator. 



the thermocline. As T approaches unity, the waves in the hurricane wake penetrate 
through the first few cycles with very little phase change in the velocity. The values of 
2.1 and 0.9 for the Norbert indicate that fairly large phase changes through the 
thermocline should be expected. The thermocline Burger number of 0.27 for Frederic is 
consistent with a deeper oceanic response that was demonstrated by Shay and Elsberry 
(1987a). The phase changes through the thermocline derived from the Frederic data 
are small with a simple vertical structure. 

Finally, the Rossby numbers, estimated from the expression T/hjU^f are 
approximately 0.2 for Frederic and 0.5 for Norbert. Chang and Anthes (1978) noted 
that the nonlinearities decrease as the translational speed of the storm increases, which 
is consistent with the scaling arguments of Greatbatch (1984). 

2. Velocity and Temperature Scales 

Price (1984) scales the maximum wind-driven velocity as xR^^/hjU^. Using 
the values in Table 4, the predicted wind-driven velocity for the hurricane Frederic case 
would be 38-40 cm/s, which is roughly one-half of the observed near-inertial velocity in 
the mixed layer (see Fig. 4.5). Similarly, the estimated mixed layer velocity during the 
passage of Norbert is about half the maximum observed velocity of 1.1 m/s (see Fig. 
5.2). The vertical velocity, which scales as t/U^, is generally larger for the slower 
moving storms because the amount of vorticity injected into the mixed layer increases. 
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For example, the predicted vertical velocity scale for Frederic is 0.06 cm/s as compared 
to 0.08 cm/s for Norbert. Furthermore, increased vertical advection in the thermocline 
displaces the isopycnals upward for the slowly moving storms. The isopycnal 
displacement T/fU^ predicted for hurricane Norbert is about 18 m, as compared to 8 m 
predicted for Frederic. Price refers to this phenomenon as "inertial pumping," which is 
a time-dependent manifestation of Ekman pumping by virtue of T/L arguments 
(dimensions of the wind stress curl). It should be noted that the time available for 
mixing also increases for slower moving storms, although mixed layer depth changes 
due to mixing are counteracted by upwelling induced by the curl of the wind stress. 

The Greatbatch (1984) scaling arguments are also used to estimate the wind- 
driven velocity, vertical velocity and temperature perturbations. As in the Price 
treatment, the Greatbatch scaling arguments apply to the area bounded by r < R x , 
and the maximum wind-driven velocities are scaled as tL/hjU^, where L is the scale of 
the wind stress curl (2R max ). Thus, the maximum horizontal velocities are 0.8 and 1 
m/s in hurricanes Frederic and Norbert, respectively. The Price and Greatbatch 
formulation of the wind-driven velocity (u') difFer by a factor of two in the numerator 
because Price uses the radius of maximum winds R max , which is half the scale of the 
wind stress curl (L) used by Greatbatch. Thus, the predicted Frederic velocity is nearly 
equal to the observed near-inertial velocity, although these data were obtained at about 
4 R max . The vertical velocities (w') are scaled similarly in both models and are 
consistent with the observations in the Frederic and Norbert storms of 0.06 and 0.08 
cm/s. 

The predicted temperature changes in the thermocline, according to the 
Greatbatch scaling, are about 0.2 °C, which is equal to the lower limit of the observed 
demodulated, near-inertial temperature changes (see Fig. 4.2) during the passage of 
hurricane Frederic. If the wind stress is nearly doubled, which is consistent with the 
Black (1983) estimates for Frederic, the temperature changes are doubled and are closer 
to agreement with the observations. The predicted temperature changes at 200 m in 
hurricane Norbert are 0.5 °C. 
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III. LINEAR THEORY 



The near-inertial velocity profiles acquired during the passage Norbert indicate 
the presence oflarge vertical wavelength oscillations. Since these profiles were 
deployed within the direct forcing regime or near-field (r < 3R max ), 11 c ^ ear t ^ iat 
vertical structure is not just a sum of the free baroclinic modes as in the far-field (Gill, 
1984). Thus, the current profiles should be examined in terms of the applied 
atmospheric forcing, especially for comparisons in the near-field. 

A. GOVERNING EQUATIONS 

The two-layer model of Geisler (1970) is extended to a continuously stratified 
ocean with a Brunt-Vaisala frequency (N(z)) following the development by Kundu and 
Thomson (1985), who treated an atmospheric front rather than a tropical cyclone. 
Geisler (1970) only derived expressions for the vertical velocity and isopycnal 
displacement. Because the AXCP's measure the baroclinic components of the 
horizontal velocities, expressions for the structure of the horizontal velocities are 
derived here. It is assumed that the ocean is linear and inviscid with a mixed layer 
depth h where N(z) = 0. The governing equations for a Boussinesq fluid are 



du 

— -fv = 

at 


1 dp 

p„ 


+X, 


(3.1) 


dv 

— +fu = 
dx 


1 dp 

p„ s y 


+ Y, 


(3.2) 


du dv 

+ + 

dx dy 


d w 

dz = °’ 




(3.3) 


dp 

dx 


Po N 2 w =o> 
g 




(3.4) 


d p 

t-+ Pg 
d z 


= 0, 




(3.5) 



where X and Y are the forcing functions (t x /hp o , t y /hp o ) distributed over the mixed 
layer as a body force. The remainder of the variables are defined as in other 
treatments. The usual kinematic boundary conditions are prescribed with zero vertical 
motion at the sea surface (z = 0) and bottom of the profile (z = -D). 
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B. FREE MODES 



Consider the free mode case in which the forcing functions are zero. Equations 
(3. 1-3.5) can then be simplified to 



The wavelengths derived from the horizontal wavenumber in the free mode case refer 
to those moving with the disturbance and would be most appropriate in the far-field. 
Since the AXCP's only sense relative or baroclinic motions, the barotropic mode can 
not be resolved from the data. The vertical structure equation is numerically solved 



order scheme (Gerald, 1983). The first five baroclinic modes for the horizontal velocity 
eigenfunctions are shown in Fig. 3.1 and the relevant parameters for the first five 
baroclinic modes are given in Table 6. 

1. WKBJ 2 Normalization 

As near-inertial waves propagate into the interior of the stratified fluid, they 
can experience variations in amplitude and vertical wavelength because of the 
stratification changes (Leaman and Sanford, 1975). The high-order baroclinic modes 
(small vertical wavelength) are more susceptible to these stratification changes. To 
accommodate these changes in the stratification, the current profiles are scaled by 



where u w is the WKBJ-scaled current at depth z w and N q is the vertically-averaged 
N(z) which ranges between 2-3 cph, so that the final scaled depth equals the original 
depth. The expressions in (3.7) represent WKBJ short-wavelength approximation to 
the wave equation (Robinson and Silvia, 1981). 

2. Baroclinic Time Scales 

If the scale of the wind stress curl (68 km) in Norbert, is greater than the 
deformation radius associated with the first mode (43 km) (large-scale limit), there will 
be a characteristic time scale for each baroclinic mode to separate from the solution in 
the mixed layer (Gill, 1984). The time required for a phase difference of 71/2 to develop 

2 The letters WKBJ stand for G. Wentzel, H. A. Kramers, L. Brillouin and H. 
Jefferys who independently discovered the procedure working on different problems 
(Morse and Fesbach, 1953). 




(3.6) 



where cp n is the n^ eigenfunction, c n = (< 7 2 - f 2 )/K n 2 is the eigenvalue and represents 
the phase speed of the n^ baroclinic mode (K n 2 = k 2 + l 2 ) and (p n = 0 at z=0,-D. 



using the observed spatially-averaged N 2 profile (Fig. 2.11) and a Runge-Kutta 4^ 



u w (z) = u(z)/V(N(z)/N o ), and 
dz w = (N(z)/N o )dz , 



(3.7) 
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Figure 3.1 Amplitudes of the horizontal velocity eigenfunctions for the first five 
baroclinic modes: 1 (solid), 2 (dashed), 3 (dotted), 4 (chain-dotted) and 5 (chain- 
dashed), based on the Brunt-Vaisala frequency profile in Fig. 2.11 and <r= 1.1 8f.. 



for each baroclinic mode n is given by 

Jtf 



<n = 



k 2 c 2 
s n 



(3.8) 



where k $ is the inverse storm scale (2R max )’^ and the other terms have been defined 
above. These time scales are given in Table 6 for the first five baroclinic modes. The 
essence of Gill's theory is that the modes will separate from the remainder of the mixed 
layer (40 m, N = 0) solution over these time scales as energy propagates vertically into 
the thermocline. This modal separation was noted in numerical simulations of the 
passage of an atmospheric front by Kundu and Thomson (1985) and for the first time 
in observations acquired during the passage of hurricane Frederic (Shay and Elsberry, 
1987a). Since the low-order modes can penetrate into the thermocline very quickly, 
these time scales will also be compared to amplitude changes in the thermocline (upper 
250 m). 
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TABLE 6 



Phase speed (c n ), deformation radius (a"^), 
wavelength (X n ), the equivalent forcing depth (j(p nz 2 (z)dz) and baroclinic 

time scales (t n ). * D 



mode 


C n 


V 1 


K 


J<p nz 2 (z)dz 


tn 


n 


(m/s) 


(km)"l 


(km) 


(m) 


(IP) 


Norbert 


1 


2.1 


43 


231 


98 


1.2 


2 


1.4 


29 


130 


334 


2.6 


3 


0.9 


19 


78 


336 


4.8 


4 


0.7 


15 


56 


421 


10.5 


5 


0.5 


11 


49 


592 


20.6 


Frederic 


1 


3.0 


42 


255 


86 


1.1 


2 


1.1 


15 


92 


174 


7.5 


3 


0.7 


10 


58 


373 


18.8 


4 


0.5 


7 


41 


651 


36.8 


5 


0.4 


6 


33 


1056 


57.5 



( ( ) nz =d( Pn/ dz 



C. FORCED MODES 

The variables in equations (3. 1-3.5) are expanded in terms of normal modes 

w,p(x,y,z,t) = X w n*Pn^ z ) ’ 

^ dtp 

u,v,p(x,y,z,t) = X u n ,v n ,p n — " . 

The variables (u,v,p,w,p) with n subscripts represent horizontal structure functions in 
(x,y,t), whereas the vertical structure functions are the q> n 's. Although the <p n 's are 
orthogonal, there is no guarantee that the d<p n 's/dz are orthogonal. The inner products 
of the d(p n /dz and the d(p m /dz are 0(0.01) whereas the inner products of (dq> n /dz ) 2 are 
0(100) as shown in Table 6. Thus, these terms are neglected in the subsequent 
analysis. Equations (3. 1-3.5) expanded into the normal modes become 



du 

— n -fv = - 

dt n 


i_ £Pp 
p„ e * 




(3.10) 


dv 

— 11 +fu = — 


J_ £En 


+ Y , 


(3.11) 


dt n 


p„ e y 


n’ 
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du dv n 

— n + _n + w = 0 , 

dx dy 



8 ±n 

dt 

Pn = 



p rt N 2 w„ 

^ n = 0j 



n 

N : 



2 §Pn 



(3.12) 

(3.13) 

(3.14) 



where the forcing functions are also expanded as ' 



T X 

X_ — 7G 



PoRz dz 



and 



Y_= 73 ” 



" PoKz 



-o 



where <p nz = d<p n /dz. 

In this treatment of the wind stress superposition onto the modes, the vertical 

c° 7 

integrals of the eigenfunction squared (J<P nz dz) will be referred to as the equivalent 

•o n 

forcing depths (Wunsch and Gill, 1976). These values are also given in Table 6. 
Notice the last two equations may be recast into 



^Pn f 1 

w = — n 

n dt \ c 



2 P 

n 



which renders the conservation of mass equation as 



du dv I dp n 

dx dy J dt 



PoV ^ ‘ 1 + T- l 7 + -e* = U • (3-15) 

Thus, the governing equation for the vertical velocity for the n^ 1 baroclinic mode 



becomes 




1 

c 



d 2 2 i 

— ^ - a 1 > w = fVxx + 



n 



2 ^ t 2 n 



d 

— V*T 
dt n 



(3.16) 



where Vxt = dY /dx - dX„/dy and V*t = dXJdx +dYJdy . The details of the 

derivations for the vertical and horizontal velocities are given in the Appendix A. 

Equation (3.16) is analogous to equation (10) in Geisler (1970) except that the phase 

speed (c n ) and the inverse deformation radius (a ) are for the n in baroclinic mode. If 

the storm is moving steadily with speed U h in the + y direction, the local time 

derivative is transformed into a space derivative 

d d 

— = -U. — . 
dt h dy 

Substituting this expression into (3.16) yields 
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where the operator is given by 



™ U. 2 d 2 d 2 

^ ( 7^2- - - 



n 



dx x 



+ a 



n 



This operator is hyperbolic (elliptical) depending on the Froude number U h /c n > ( < ) 
1. Only the hyperbolic case (U h > c n ) is considered here. Since the AXCP's measure 
the vertical shear of the horizontal velocity, the spatial distribution of the 
corresponding eigenfunctions (u n (x,y), v n (x,y)) is sought using (3.10), (3.11) and (3.15), 
which yields (Longuet-Higgins, 1965) 



<2> n = - + -L (fY - 



dX n 

2 " n_ ~~ U^ VA '' n ' V 2 Vlln Lfl dv ^ 

n 



'h 

_ \ f d 1 dY n 

= — f — Vxrdy 5 (fX+U. — n ), 

2 n U h - J oo dx n - 2 n h '• 



n 



dy 



where the integrand in the v n equation is defined only in the direct forcing region. The 
details of this derivation are given in Appendix A. Geisler (1970) introduced a 
nondimensional coordinate system 



y = 



Hu 

c „ 2 



-1, - a 



■ n y» 



x - a n x 



which transforms the operator into 

V <3 2 



££ Y = a 2 

Z n 



dy 2 dx ’ 2 + 1 



The governing expressions for the vertical and horizontal velocities are 



™ 1 U h d 

££Y\v = Vxt + — 

2 n f n 



K 3T' 7,t " ) ’ 



dy 

1 1 dX n 

«£?,u = Vxt + ^ (fY-U. — n ), 

2 n U. a 2 n f 2 n h dv 

n n 



(3.17) 



2 n 



1 ~d fc 7xT " dy " ? (K " +U *> 



U 



h n 



d Y_ 
dy 






The VxT n and V*T n are the wind stress curl and divergence of the wind stress as a 
function of mode number. The dominant wind stress curl term sets up the semi- 
permanent baroclinic ridge in the wake of the storm. The second term in the vertical 



velocity equation is the gradient of the wind stress divergence. Since the divergence 
term is negative (convergence), the gradient of the divergence has the same sign as the 
wind stress curl, which augments the effect of the wind stress curl on the ocean current 
response. In the horizontal velocity equations, the second terms on the right side 
represent those components that are directly driven by the wind stress (Y /f and -X n /f), 
which are referred to here as Ekman-type velocities. The third terms in the u and v 
equations represent a sink of vorticity and the y-component of the divergence, 
respectively. Most analytical treatments have ignored the divergent and Ekman-type 
terms because the wind stress curl sets up the semi-permanent baroclinic ridge, and it 
has been assumed that the inflow angle was zero. However, Mayer et al. (1981) 
showed that the scales of the wind stress divergence were 30 - 50 % of the curl term in 
hurricane Belle. 



D. IMPULSIVE FORCING 



Since the translation speed of Norbert exceeds the first mode internal wave phase 
speed, the oceanic response should be dominated by baroclinic near-inertial waves in 
the wake of the hurricane. A further consequence of this relationship is that the ocean 
is undisturbed well ahead of the storm, which may be demonstrated by a direct 
application of the radiation boundary condition (Lighthill, 1967). 

The solutions of (3.17) begin with the assumption that the right side can be 
approximated with a dirac delta function to represent an impulsive type of forcing 
function 

Vxt n = -t on L 5(x)5(y) , 

where T Qn represents the magnitude of the wind stress projected onto the n 1 * 1 baroclinic 
mode and L is the the scale of the storm (typically 2R max ). Without loss of generality, 
only the curl term will be operated upon here, and the divergence and Ekman-type 
terms will be dropped from the transform and inversion process (but not from the 
solution). The implusive forcing assumption can be similarly applied to the other 
terms. 



The first equation in (3.17) is recast into an expression with an impulsive forcing 
on the right-hand side with the nondimensional operator on the left side 
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Double Fourier transforming (3.17) yields 
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and substituting into the inversion formula renders an integral equation representing a 
fundamental solution to the problem or a Green's function (Stakold, 1978) 



t a ^ r P ^ i oo 

F (x',y) = -f- 1 A -i} 1/2 ft 

tU C J .oo 



1 



,-ikx' e -ily' dk dl _ 



(k 2 -l 2 + 1) 

The Green's function that emerges from the solution of the integral equation is a 
Bessel's function of order 0 (see Appendix B): 

1 ikx' ilv' J n ((y' 2 -X' 2 ) 1/2 ) 

f lkx e ll y dk dl = -2^ — - 



00 

ft 

•00 



(k 2 -l 2 + 1) 



H(y'-x') . 



Because a hurricane can be represented by a circularly symmetric point source 
translating over the ocean, the Bessel function is geometrically appropriate and H(y'-x') 
represents a unit step function. Since the argument of the Bessel function involves 
(y',x') or (a y, a n x), the fundamental solution also represents a two-dimensional 
eigenvalue problem. 



E. THE SOLUTION 

Generally, the solution to an integral equation is represented by a convolution of 
the forcing term and the Green's function integrated over all source points (T,r) 
within the direct forcing region 

«'„= - jj 7xt n (r,y)G(y',r;x,r) dr djr , 

X «oo 

on 

where the coefficient (T on L)*^ is due to the impulsive forcing assumption. 

As noted above, the complete set of forcing functions will now be introduced into 
the convolution integral for the 3-dimensional velocity structure 



w „- v Ii~ x W r ') d *' d r - $r 

11 • oo O ll-co O o y 
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00 ^* 00 - 00 o cy 

where the Bessel function represents the kernel in the convolution and the argument is 
given by: 
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The limits of the integration are from zero to -oo in y' as a result of the radiation 
boundary condition. In these convolutions, the direct forcing region will be restricted 
to R max for the curl and divergence terms, and to 40 km for the Ekman-type terms. If 
the size of these regions is doubled, the simulated velocities are similarly increased. 



F. COMPARISON TO THE TWO-LAYER MODEL 



The effect of a continuously stratified ocean model versus a two-layer model as in 
Geisler (1970) is demonstrated with a wind stress curl similar to that of Geisler 



Vxr n = 2.5 



ijl- .iU 



(3.20) 



where T m is the maximum wind stress (3 N/m 2 ), r is radial distance from the storm 
center and L is the length scale of the wind stress curl (100 km). In these numerical 
experiments, the storm translation speed is 7 m/s and a phase speed of the internal 
mode is 1.4 m/s. 

In this version of the model (Fig. 3.2a), maximum upward velocities are slightly 
greater than the 3 nondimensional units predicted by Geisler (1970). In dimensional 
units, these values are about 7-8 x 10'^ cm/s, which is twice the velocity Price (1983) 
and Greatbatch (1984) predicted from scaling arguments. The pattern of vertical 
velocities is elongated in the cross-track direction over the first half of the inertial cycle 
(220 km) with a scale of 8R max . The lateral spreading of the wave-wake is due to the 
localized core of positive vorticity, which is usually confined to r < 2R max . 

To compare the effect of various wind stress curl patterns, a second numerical 
experiment is performed by replacing the stress pattern in (3.20) with a Rankine vortex 
(2.2) that is based on the same parameters above. The region of positive vertical 
velocities (Fig. 3.2b) is confined to a narrower region of r < 2R max with comparable 
maxima of 7-8 x 10 . The pattern is more crescent-shaped and more resembles the 
character of the Bessel function than in the simulation with (3.20). The vertical 
velocities damp more quickly because the spatial scale of the wind stress forcing is 
smaller by 2R max . 

1. Vertical Velocity 

Using the Rankine vortex and storm parameters from Norbert, the total 
vertical velocity for a sum of the first three baroclinic modes is shown in Fig. 3.3. The 
vertical velocity is upward (positive) in the direct forcing regime with maximum values 
of 0. 1-0.2 cm/s which agrees well with scaled vertical velocity (0.1 cm/s) based on the 
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maximum stress. The upward (downward) vertical velocities indicate areas of 
upwelling (downwelling) on the semi-permanent baroclinic ridge that is established due 
to the net input of cyclonic vorticity by the storm (Price, 1981). Upper ocean currents 
diverge from the storm track as cooler water is upwelled. Conversely, the horizontal 
currents in areas of downwelling, converge toward the storm track. The horizontal 
wavelength for the first mode estimated from linear theory is 440 km. This is the 
wavelength in the fixed frame as opposed to the wavelength moving with the storm, 
w'hich scales with wavelength of the first free mode^Cj = cr/Kj). 

As expected, the pattern is dominated by the wind stress curl (Fig. 3.3a). For 
an inflow angle 20 °, the contribution of the divergence term in the direct forcing 
region is approximately 25 % of the curl term. For larger inflow angles, the w r ind 
stress divergence becomes larger. However, there does not seem to be any semi- 
permanent ocean features set up by the wind stress divergence as in the case of the 
curl. The effect of the divergence damps away from the directly forced region within 
2-3 inertial wavelengths. The contribution here is associated with the gradient of the 
divergence, which has a scale half that of the curl term and accounts for the decrease in 
magnitude. 

Since the kernel of the integral in (3.19) is the Bessel function, the decay of the 
velocity amplitudes follows from an asymptotic analysis of that function (Arfken, 1970) 

J Q (y)= V— cos (y- 71/4) 7 (3.21) 

which implies that the amplitude approximately decreases as y'^. 

The structure of the horizontal velocities (Fig. 3.4) for a summation of the 
first three modes corresponds to the pattern of the vertical velocity. Maximum v- 
components exceed 100 cm/s in the first upwelling area in the wake. Within the next 
half-wavelength, the v-component changes direction and is associated wdth the 
convergence of flow and the downwelling of the near-surface waters. The u-component 
velocity pattern is out of phase with the v-component and the maximum u-component 
is 60 cm/s. Thus, the velocities are varying as expected for forced, near-inertial 
motions superposed on the mean flow associated with the baroclinic ridge. The pattern 
broadens in the wake of the storm according to tan'* (TJ h 2 /c n 2 -l) - ^ ^ and the 
amplitudes decrease according to the asymptotic expression (Geisler, 1970). 

The major term in the momentum balance of the horizontal velocites is the 
wind stress curl (Fig. 3.4b). As expected, the first baroclinic mode terms dominate the 
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Figure 3.3 Nondimensional vertical velocity contours at 90 m for a summation of the 
first three baroclinic modes in the thermocline based on the Rankine vortex model for 
(a) only the wind stress curl, (b) only the wind stress divergence and (c) the total 
(a + b). The model storm is moving in the 4-y direction at 4 m/s. The contour interval 
is as Fig. 3.2 and the ordinate and abscissa axes arc nondimensionally scaled with the 
wavelength of the first mode (440 km) and R m (30 km). 
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other baroclinic modes by factor of about three to four because of the superposition of 
the wind stress onto the modes. The divergence and Ekman terms are considerably 
weaker with relative maxima of 20 and 10 cm/s, respectively (Fig. 3.4c and d). The 
divergence term augments the wind stress curl term in the direct forcing area, but 
damps very quickly after 2 A. By contrast, the Ekman terms act to diminish the 
velocities induced by the curl and divergence terms by roughly 10 cm/s. The 
importance of the Ekman term also diminishes rather quickly, w r hich is not surprising 
since the wind stress is also diminishing quickly beyond a radius of 200 km. It is clear 
that the wind stress curl dominates the ocean response over short and long space (time) 
scales. However, the divergence and Ekman terms do contribute to oceanic current 
variations within the direct forcing region. 

2. Simulated Vertical Structure 

Using the expansions in (3.9), the simulated velocity fields from the linear 
model (u n (x',y'),v n (x',y')) are constructed in the form of a depth-space or depth-time 
series by multiplying the horizontal structure functions by the vertical structure 
functions (p(z). Since the model is linear and the large scale limit is satisfied (S= 1.5), 
these separation time scales of the modes predicted by Gill (1984) should become clear 
in the wake. As shown above (Table 6), these baroclinic time scales are : t 1 = 1.2 IP; 

U = 2.6 IP; and t 3 = 4.8 IP in hurricane Norbert. The essence of Gill's theory is that 
the modes will separate from the remainder of the mixed layer solution over these time 
scales as energy propagates vertically into the thermocline. This modal separation was 
noted in numerical simulations of the passage of an atmospheric front by Kundu and 
Thomson (1985) and for the first time in observations acquired during the passage of 
hurricane Frederic (Shay and Elsberry, 1987a). 

An along-track (or time series) section of the v-component for a summation of 
the first three baroclinic modes is shown in Fig. 3.5 at x= ± R max - The abscissa 
depicts either space or time normalized to a wavelength A of 440 km or an IP of 1.5 
days. Within the first 1.5 IP, the maximum velocity component decreases from about 
110 cm/s to about 70 cm/s as the first baroclinic mode propagates out of the mixed 
layer and into the thermocline. The slope of the zero lines is upw r ard towards the left 
which represents an upward phase propagation. This corresponds to downward energy 
propagation for near-inertial waves (Leaman and Sanford, 1975). A second decrease in 
mixed layer velocity after about 3 IP corresponds to the time scale associated with the 
second baroclinic mode. Energy in the upper ocean remains essentially constant 
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Figure 3.4 Simulations of the mixed layer v-component (cm/s) for a summation of the 
first three baroclinic modes for forcing due to (a) all the forcing terms, b) the wind 
stress curl, (c) the wind stress divergence and (d) the wind stress based on the observed 
storm parameters in Norbcrt. The contour interval is 10 cm/s with a maximum of 110 



cm/s. 
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between 3-5 IP. The time required for the first baroclinic mode to again become in 
phase in the mixed layer is roughly 4tj or 4.8 IP, which also corresponds to the time 
scale of the third baroclinic mode separating from the mixed layer (t 3 = 4tj). At about 
5 IP, the mixed layer velocity actually increases by 10 cm/s, which appears to be 
associated with the more dominant first baroclinic mode rather than a decrease in the 
mixed layer velocity as the third mode separates from the solution. This full-cycle of 
mode 1 (4tj) also agrees with simulations by Kundu and Thomson (1985). 

G. SUMMARY 

In this chapter, the linear, two-layer model of Geisler (1970) has been extended 
to include a continuously stratified fluid. The wind stress is superposed unto the 
first five baroclinic modes as in Kundu and Thomson (1985). The veracity of the 
analytical model was demonstrated by collapsing the continuously stratified model into 
a two-layer model to directly compare with the results of Geisler. The slight differences 
in the response of the vertical velocity is due to the selection of the forcing functions. 

It appears that the Rankine vortex forcing gives more realistic results than the forcing 
(3.20) used by Geisler. 

To facilitate comparisons to the AXCP observations, explicit relationships for the 
horizontal velocities have been derived for a northward moving storm. The behavior of 
the modal amplitudes for the v-component agrees well with the time scales predicted 
from Gill's theory. Thus, the analytical model velocities for the first few modes "will be 
compared to the AXCP data acquired during Norbert within the near-field. 
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Figure 3.5 Summation of the first three baroclinic modes for the v-component from 
the forced model in the along-track direction at x=R max (upper) and 
x* -R max (lower) with a contour interval of 10 cm/s. The solid (dashed) contours 
represent positive (negative) velocity. The abscissa is normalized by the wavelength of 
the first mode (440 m). 
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IV. TIME EVOLUTION 



The advantage of the moored current meter data is the ability to depict the 
detailed time evolution of the currents and temperatures at selected levels prior, during 
and subsequent to storm passage. In this chapter, the near-inertial response is isolated 
from the moored ocean current and temperature measurements in hurricane Frederic 
by complex demodulation. The amplitudes of the near-inertial currents are fit to the 
the vertical modes and the time-evolution of the modes is compared to the predictions 
from linear theory of Gill (1984). Because these data were acquired along the 
periphery of the DeSoto Canyon, the effects of the sloping bottom on the first two 
baroclinic modes are studied using the model of Lai and Sanford (1986). 

The ocean current and temperature data were low-pass filtered at 3 h using a 
Lanczos taper window (Table 7) and were subsampled every hour to smooth the data 
series. The smoothed data series were used for most of the data analyses except in the 
spectral calculations where increased time resolution was required. 





TABLE 7 


Characteristics of the low-pass Lanczos filter. 


Period (h) 


Energy Rejection 


3 


6 db at 2.9 h 




20 db at 2.6 h 


48 


6 db at 40 h 




20 db at 24 h 



The /nertial Periods (IP's) are equal to 24.2 and 24.5 h at CMA 2 and 3 (Fig. 

2.1), respectively whereas the local IP at the OTEC site is 25.8 h. The length scale of 
the near-inertial response will be referenced to the scale of the wind stress curl, which is 
approximately 60 km for Frederic (Powell, 1982), and the Rossby radii of deformation 
(see Table 6) for the barotropic and first three baroclinic modes. These deformation 
radii are based on the near-inertial wave phase speeds computed from the Sturm- 
Liouville problem, which was discussed in Chapter III. 
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If the scale of the wind stress curl (60 km) is greater than the deformation radius 
associated with the first mode (42 km), there will be a characteristic time scale for each 
baroclinic mode to separate from the solution in the mixed layer (Gill, 1984). For the 
Frederic case, the baroclinic time scales are tj = 25 h, t2 = 7.5 IP and t^ = 18.8 IP 
for first three baroclinic modes. The above relationship was derived for changes in the 
baroclinic modes within the mixed layer. Since the low-order modes can penetrate into 
the thermocline very quickly, these time scales will also be compared to clockwise 
amplitude changes in the thermocline (upper 250 m). 

A. OCEAN CURRENT SPECTRA 

A Tukey data window (Otnes and Enochson, 1978) is applied to the original 
current observations prior to the removal of the mean in the calculation of the energy 
spectra. The data are Fourier transformed and spectrally-averaged over bandwidths to 
help decrease leakage of energy to adjacent frequency bands (Otnes and Enochson, 
1978). Energy above the Nyquist frequency (l/2At), where At is the 10 minute 
sampling interval ( 20 minutes at OTEC), is eliminated to minimize aliasing of the 
spectra. 

The auto-spectra of the cross- and along-track velocity components in the mixed 
layer at CMA3 are fairly representative of the spectra throughout the water column 
(Fig. 4.1). The dominant peak in the spectra is shifted above the local inertial 
frequency and is significant at the 95 % confidence level. The variance associated with 
the near-inertial motion is nearly two orders of magnitude larger than the peak at the 
semi-diurnal frequency. Because of the large bandwidth in the near-inertial band 
(typically 0.01 cph), the frequency shifts cannot be determined from the spectra. The 
spectra are red at the higher frequencies, which indicates fairly reliable current 
measurements. 

The energy spectra are decomposed into clockwise (CW) and counter-clockwise 
(CCW) rotating components to illustrate the rotational characteristics of the inertial 
wave motion (Gonella, 1972 ; Mooers, 1973). The rotary spectrum estimates for the 
near-inertial period motions are given in Table 8. Generally, the CW spectral density 
estimates are 1 to 2 orders of magnitude larger than the CCW rotating motion. 
However, the CCW spectral density exceeds that in the CW direction in the upper 
thermocline (100 m) at the OTEC site. Since, the CW and CCW spectral densities at 
CMA1 are not significantly different, the paper will concentrate on the near-inertial 
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Figure 4.1 Normalized auto-spectra from CM A3 for a) cross-track and b) along-track 
velocity components in the mixed layer. Kj, f, M2 represent the diurnal tide, inertial 

and semi-diurnal tide Irequencies. The 95 % confidence intervals are shown at the 
bottom. 



response at CMA2, CMA3 and the OTEC moorings. The rotary coefficients are 
indicative of the amount of coherence or correlation between the horizontal velocity 
components (Fofonoff, 1969). For near-inertial oscillations, the coefficients approach 
unity, which is the case in these records except in the upper thermocline at OTEC. By 
contrast, the subinertial (periods of 3 IP) and superinertial (periods of 12 h) CW 
spectral estimates (not shown) are more than an order of magnitude less than the near- 
inertial motion. 

1. Tidal Analysis 

A harmonic analysis was performed over a 29-day time series of the pre- 
hurricane current data following the method of Doodson and Warburg (1941). The 
variances associated with the semi-diurnal and diurnal frequency bands are given in 
Table 9 for the pre-storm period. Generally, the tides in this region are diurnal with 
very weak semi-diurnal components (Zetler and Hanson, 1972). 

The diurnal tides are dominated by the Kl, 01 and PI tides in this region of 
the Gulf of Mexico (Molinari et al . , 1979). Considerably more variance is associated 
with the diurnal tides in the upper layers than in the lower layers. Molinari and Mayer 
(1982) also found larger tidal currents at 100 m (4-5 cm/s) than at 550 m (2-3 cm/s) at 
the Tampa OTEC site. The corresponding amplitudes in the mixed layer at CMA2 and 
CMA3 were 4-5 cm/s. Below the mixed layer, the tidal amplitudes were 2 - 3 cm/s, 
except at 437 m (CMA3) where the amplitudes nearly equaled that observed in the 
mixed layer. These estimates agree with other estimates in the Gulf of Mexico 
(Brooks, 1983 ; Daddio et al., 1977). The periods of near-inertial and diurnal tidal 
signals are very close and cannot be resolved since the spectral bandwidth around the 
daily oscillations is about 0.01 cph. Thus, the diurnal tides were not removed from the 
time series. 

2. Current and Temperature Response 

Near-inertial wave excitation by hurricane Frederic was clearly evident 
throughout the water column at the three arrays (Fig. 4.2). The increase in the 
currents at all levels was first felt about 8-9 h in front of the storm. The mixed layer 
currents oscillated with magnitudes of 80-90 cm/s, which agrees with the mixed layer 
currents in the numerical simulations of hurricane Frederic by Hopkins (1982) and with 
predictions based on the scaling arguments of Greatbatch (1984). These currents then 
decreased over c-folding scales of about 4 IP. However, the behavior of the near- 
inertial waves in the thermocline was markedly different from the mixed layer. For 
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TABLE 8 



Near-inertial rotary spectrum estimates shifted above the 
local inertial frequency. 



Rotary 





Depth 

(m) 


CW 

(cm/s) 2 /cph 


ccw 


Coeff. 


CMA1 




49 


6.7x10-* 


1.3x10-* 


+ 0.68 




64 


2.5xl0 3 


1.4xl0 3 


+ 0.62 


CMA2 




19 


l.SxlO^ 


3.5x10-* 


+ 0.96 




179 


4.6xl0 4 


3.7xl0 2 


+ 0.95 




324 


2.8xl0 4 


1.3xl0 2 


+ 0.96 


CMA3 




21 


4.0xl0' ) 


2.0x10-* 


+ 0.98 




251 


1.4xl0 5 


6.1xl0 2 


+ 0.98 




437 


1.2xl0 5 


6.3xl0 2 


+ 0.98 




457 


6.2xl0 4 


5.0xl0 2 


+ 1.00 


CMA6 




19 


lAxlO* 


5.1x10-* 


+ 0.95 




180 


7.5xl0 3 


2.0xl0 2 


+ 0.92 




330 


5.0xl0 3 


l.OxlO 2 


+ 0.95 


CMA9 




20 


3.8x10^ 


1.8x10-* 


+ 0.95 




250 


l.lxlO 5 


5.0xl0 2 


+ 0.93 




500 


5.0xl0 4 


2.0xl0 2 


+ 0.95 


OTEC 




100 


2.7x10-* 


3.3x10-* 


-0.20 




232 


2.1xl0 4 


l.OxlO 2 


+ 0.95 




546 


1.6xl0 4 


9-OxlO 1 

7.0x10* 


+ 0.95 




950 


3.5xl0 3 


+ 0.90 



example, the currents in the thermocline at CM A3 (251 m) and CMA2 (179 m) 
oscillated with amplitudes of 25 cm/s and described a modulation envelope within 7 IP 
subsequent to the passage of Frederic. The thermocline currents at CMA3 persisted 
for nearly the entire record following storm passage (21 IP). The near-bottom currents 
showed similar modulation behavior at both arrays. The OTEC data (Fig. 4.2c) 
indicate that the near-inertial response is even deeper than suggested by the 
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TABLE 9 



Variance of the semi-diurnal (M2,K2,N2,T2) and diurnal 
(K 1,01, PI) tides estimated from 29 days of pre-storm data. 





Depth 


Semi-diurnal 


Diurnal 




(m) 


(cm/s)^ 


(cm/s)^ 


CMA1 




21 


3 


44 




49 


2 


5 




64 


2 


6 


CMA2 




19 


2 


23 




179 


2 


10 




324 


1 


1 


CMA3 




21 


2 


20 




751 


1 


4 




4 37 


1 


16 




457 


1 


3 


OTEC 




100 


2 


# 




232 


1 






546 


1 


M. 




950 


2 


s 

7T 



* OTEC analvses are based on only 10 days of pre-storm data. 

# Could not 6e resolved from 10 days of pre-storm data. 



NAVOCEANO data. Although the current speeds at 232 m were considerably weaker 
(±15 cm/s) than those observed at the CMA3 (which was closer to the hurricane 
path), the currents at the OTEC site increased in magnitude between 4-6 and 8-10 IP. 

The thermocline and near-bottom temperatures at CMA2 and CMA3 were 
modified considerably by the passage of the storm. These temperatures were in-phase 
with the along-track velocity component, as expected from plane wave theory (Mayer 
et al., 1981). At CMA2 (Fig. 4.2b), the pre-storm temperatures at 179 m ranged 
between 14.3 to 14.8 °C. As the storm arrived, the temperature increased by about 
3 0 at this level. This was followed by a near- inertial response over the next 10 IP. 
Initially, the near-bottom (324 m) temperature decreased by 1.5 0 C and subsequently 
increased about 3.5 0 C over one IP. Small amplitude, near-inertial oscillations were 
superposed on a lower frequency (about 4 IP) oscillation following storm passage. At 
CMA3, the thermocline (250 m) and near-bottom (457 m) temperatures followed 
similar trends except for a shorter period oscillation (2 IP) in the near- bottom 
temperatures. 
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The near-inertial response in the ocean currents and temperatures at CMA6 
and CMA9 are quite similar to those at CMA2 and CMA3. The thermocline currents 
at CMA9 (250 m) increased in amplitude to about 22 cm/s after about 7 IP which also 
defined a modulation envelope observed at CMA2 and CMA3. 

Mean currents were estimated by a running average over an inertial period. 
The largest mean currents were oriented along the bottom topography at all three sites. 
At CMA3 (Fig. 4.3a), the mean along-track (nearly parallel to the topography) current 
direction was towards the south at 4 cm/s. However, there was a significant change in 
the mean current about 9-12 h prior to the time of closest approach of hurricane 
Frederic. First, the mean currents increased at all four depths by about 10 cm/s 
towards the north. Subsequently, the mean current reversed over 2 IP and exceeded 10 
cm/s towards the south. Moreover, the near-bottom mean current exceeded that in the 
mixed layer and thermocline. The mean cross-track current in the mixed layer at 
CMA2 was 20-25 cm/s before the arrival of the storm and was much stronger than at 
CMA3 (Fig. 4.3b). Although the mean current decreased slightly during the period of 
inertial wave excitation, it returned to pre-storm values within an IP. The mean mixed 
layer current reversed after 8 IP, but it is unclear whether this was related to the 
forcing associated with hurricane Frederic. Another significant feature at CMA2 is the 
increase of the thermocline mean current to 25 cm/s about 5 IP following the storm. 
Although the mean current at this level did not reverse as in the mixed layer, a rapid 
decrease to pre-storm values occurred after 8 IP. The observed mean currents at 
CMA6 and CMA9 follow the bottom topography and support these interpretations. 

The mean flow at the OTEC site (Fig. 4.3c) also changed a few hours prior to 
the time of closest approach of hurricane Frederic (12 UTC 12 September). The mean 
currents exceeded 5 cm/s flowing towards the north at all depths. The mean current at 
100 m oscillated with periods on the order of 5 IP after storm passage. This appears to 
be a local phenomena as it was not observed at any of the other depths. 

The mean currents in the thermocline were also markedly different at all array 
sites (see Fig. 4.3). Whereas a slow, southward mean flow was observed at CMA3, the 
mean flow was towards the north at all levels at the OTEC site. This pattern suggests 
a convergence zone existed between the DeSoto Canyon region and the OTEC site. 
Climatologically, there is a predominant northward mean flow following the isobaths 
during the summer months (Molinari and Mayer, 1982). Furthermore, Black (1983) 
indicates that the Loop Current was located about 100 km south ot the OTEC site 
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Figure 4.3 Running inertial period mean of the velocity component (cm/s) following 
the bottom topography at a) CMA3 (along-track), b) CMA2 (cross-track), and c) 
OTEC (along-track) with same time scale as in Fig. 4.2. 




during the passage of Frederic. Thus, the mean circulation in this region may be quite 
complicated. 

3. Frequency Analysis 

Near-inertial oscillations are defined to be in the frequency range of 0.90f to 
1.20f where fis the local Coriolis parameter (Mooers, 1975 ; Kunze, 1985). The 
Frederic data were least-squares fit to the model of Mayer et al. (1981) to diagnose the 
period of the near-inertial oscillations. For a series of trial frequencies, a set of weights 
(uj,^) was determined from the expression 

u(t) = UjCos(at) + U2sin(<Tt) + Uj^ , (4.1) 

where <7 is the trial frequency (0.90f- 1.200 and is the residual current after the 
removal of the signal with that frequency. A similar expression was used for the along- 
track velocity component. The product of the residual currents (Ur,Vr) were averaged 
in time to form a covariance between the two velocity components at zero lag. The 
carrier frequency is defined as the frequency that minimizes the covariance of the 
residual currents. Because of the strength of the near-inertial response to the 
hurricane, it is assumed a single carrier frequency exists. The fit was performed on the 
data starting 12 h following the point of closest approach of Frederic to avoid 
contamination of the signal by the initial transient circulation (Mayer et al., 1981). 

The current data were fit to the frequencies over segments of 1, 3 and 7 IP for each 
time series. More stable estimates of the near-inertial frequency were obtained with the 
7 IP averages. 

The variance reduction curves for CMA3 are shown in Fig. 4.4. The carrier 
frequency may not necessarily be aligned with these maxima, since it is defined to be 
the frequency that minimizes the covariance of the residual signals. Near the bottom, 
the diagnosed frequencies are 1.06f (based on the minimum covariance), but the 
relative maxima in the explained variance occurs at a = 1.08f (Fig. 4.4c and d). Thus, 
there is uncertainty of about 0.02f or about 0.4 h. In the mixed layer and thermocline, 
the explained variances exceed 70 % with very little uncertainty in the diagnosed 
frequency. 

The diagnosed frequencies are generally blue-shifted between 1-6% above the 
local inertial frequency (Table 10) w’hich is predicted from the mixed layer Burger 
number (see Table 5). This blue-shift in the frequency indicates that the forced near- 
inertial waves rotate faster than the local inertial period. For these values, M is 
roughly 0.08 in the Frederic case or a predicted blue-shift in the frequency of 0.04 in 
the mixed layer, which nearly equals the diagnosed frequency at CMA2 and CMA3. 
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Variance (X) Variance (X) 





Figure 4.4 Explained variance (%) for the cross-track (solid) and along -track (dashed) 
velocity components at trial frequencies normalized by the local inertial frequency (<J7f) 
from CMA3 for a) 21 m, b) 251 m, c) 437 m, and d) 457 m. These estimates are based 
on least-squares fit over 7 IP. The vertical dashed line represents the frequencv at 
which the covariance of the residual is minimized. 
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TABLE 10 



Near-inertial frequencies (periods) at which a maximum 
response occurs in the observed currents based on a least-squares fit over 7 







IP. 




CMA2 


Depth 

( mf 


<r/f 


Period 

(h) 




19 


1.04 


23.3 




179 


1.02 


23.7 




324 


1.05 


23.1 


CMA3 


21 


1.04 


23.5 




251 


1.02 


23.9 




437 


1.06 


23.1 




457 


1.06 


23.1 


OTEC 


100 


1.04 


24.8 




232 


1.01 


25.5 




546 


1.03 


25,0 




950 


1.04 


24.8 



Since the lower order modes propagate rapidly into the interior of the fluid, a 
consistent blue-shift is expected with depth (Price, 1983 ; Kundu and Thomson, 1985). 
However, the blue-shifts in the frequencies in this case do not always increase with 
depth as the frequencies in the thermocline (typically 1.02f) are less than in the mixed 
layer. Another exception is at the OTEC site where the frequency at 232 m is less than 
at 100 m. The uncertainty of 0.02f in determining the exact frequencies as mentioned 
above can shift the frequency to that observed in the mixed layer. One possible 
mechanism for shifting the frequency of near- inertial waves is the mean flow (Mooers, 
1975 ; Kunze, 1985). If the scale of the mean flow is comparable to the wavelength of 
the propagating near-inertial waves, the relative vorticity of the mean flow' can 
effectively shift the frequency of the near-inertial waves. Thus, the frequency can be 
shifted above (below) fin regions of cyclonic (anticyclonic) vorticity. For example, 
Gaul (1967) found an anticyclonic rotating gyre in the DeSoto Canyon region that 
might account for the shift in frequency of the waves towards f. Unfortunately, the 
relative vorticity field can not be resolved with the available data in this case. 
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B. HORIZONTAL STRUCTURE 

The near-inertial component in the time series was isolated by demodulating the 
ocean current and temperature data (Perkins, 1970; Otnes and Enochson, 1978). The 
data were multiplied by the function exp (j27taiAt) , where j is +J-\, <7 is the carrier 
frequency (see Table 10) and iAt represents the time at the i^ data point. The data 
were low-pass filtered at 48 h and combined to form the amplitude and phase time 
series of the CW and CCW-rotating components (Mooers, 1973). 

1. Mixed Layer 

The increase in the CW-amplitude following hurricane passage was similar at 
both CMA2 and CMA3 (Fig. 4.5). The maximum CW-amplitude persisted for about 2 
IP at CMA3, but at CMA2 the near-inertial amplitude started to decrease immediately. 
After 4-4.5 IP, the near-inertial amplitudes at CMA2 began to increase again in the 
mixed layer and nearly equaled the near-inertial amplitude at CMA3 after 5 IP. This 
time scale roughly corresponds to the time (4tj) required for the first baroclinic mode 
to become in-phase with the higher order modes in the mixed layer (see Fig. 10c in 
Gill, 1984). This provides further observational evidence of the interchange of energy 
between the mixed layer and the stratified fluid below (Gill, 1984 ; Kundu and 
Thomson, 1985). Differences in the near-inertial response between CMA2 and CMA3 
may have been due to local mean currents. The mean current direction followed the 
bottom topography and w r as stronger at CMA2 than at CM A3 (see Fig. 4.3). The 
maximum amplitudes are about 70 and 90 cm/s after the removal of the surface wave 
current components at CMA6 and CMA9, respectively. At both arrays, a slight 
secondary increase in the CW-amplitudes is evident between 4-5 IP as observed at 
CMA2 and CMA3. 

The energy flux from the mixed layer can be estimated using the relationship 
pu^h/2T where p is the density, u is the amplitude of the near-inertial current and T is 
e-fold time scale of energy (5 IP). The energy flux estimate from the mixed layer based 
on the Frederic observations is 23 mW/m^. Brooks (1983) reported an estimate of 10 
mW/m^ in the mid-thermocline during the passage of hurricane Allen. 

Initially, the near-inertial wavelength of the ocean response is set by the scale 
of the wind stress curl in a hurricane (Price, 1983). The complex coefficients of the 
CW-rotating current vectors can be used to estimate the correlation coefficient and 
phase angle between two instruments over the first 7 IP (Kundu and Thomson, 1985). 
The phase angle is measured counter-clockwise from the CW-rotating vector at CMA2 
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100-1 




Figure 4.5 Clockwise (CW)-amplitudes (cm/s) of the near-inertial component in the 
mixed layer at CMA2 (solid) and CMA3 (dashed). Time scale is the same as in Fig. 
4.2. 

/ 

relative to that at CMA3 and is weighted by the magnitude of the instantaneous CW- 
rotating vectors. Assuming a plane wave relationship in a mixed layer of constant 
depth, the phase difference can be used to estimate spatial wavelengths (X f ) from 
dQ;dr = 2 k/X t , or 

X f = 27tAr/A0, (4.2) 

where Ar is the horizontal (vertical) distance between two instruments and A0 is the 
phase difference in radians. The spatial separation between CMA2 and CMA3 is 39 
km. 

The phase angle between the mixed layer currents is about 55-60 ^ during and 
subsequent to the passage of the storm (Fig. 4.6) and remains phase-locked during the 
period of large CW-amplitudes. The 1-h difference in time between the points of 
closest approach of Frederic to CMA2 and CMA3 has been incorporated into the 
phase. The estimated wavelength from (4.2) is 240-260 km, which is roughly 8 R max 
or four times the scale of the wind stress curl (60 km). The 250 km wavelength agrees 
with the estimated wavelength of the first baroclinic mode that is derived from the 
Sturm- Liouville problem using a variable profile (see Table 6). Black (1983) also 
observed wavelengths of 240 km in the AXBT data acquired during the passage of 
hurricane Frederic. However, linear theory by Geisler (1970) suggests that the 
horizontal wavelengths should be much larger. For a storm speed of 6.5 m/s and IP = 

24.4 h, the predicted wavelength should be about 550 km. The horizontal wavelengths 
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of the ocean response to hurricane Allen (370 km) were also less than predicted by 
Geisler's linear theory (Brooks, 1983). Because the mechanism that forces these shorter 
horizontal wavelengths of the near-inertial response is related to the wind stress curl 
(Price, 1983; Greatbatch, 1984), the wavelength associated with the near-inertial 
response varies with each hurricane. 
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Figure 4.6 CW-phase in the mixed layer currents between CMA2 and CMA3. The 
start time is 21 UTC 11 Sept, and ends 21 UTC 19 Sept, and represents the period of 
significant CW-amplitudes in Fig. 4.5. 

The CW-rotating phases (Table 1 1) are compared at the various array sites. 
These additional arrays provide five more estimates of the horizontal wavelengths of 
the near-inertial response, thus increasing the confidence in the estimate of the near- 
inertial wavelength based on CMA2 and CMA3. In these other five cases, there is a 
phase-locking of large CW-rotating amplitudes over the first 5 IP in the mixed layer (as 
in Figs. 4.5 and 4.6). The corresponding estimates of wavelengths are comparable to 
the 250 km wavelengths estmated between CMA2 and CM A3 and range between 
200-300 km. Thus, the wavelength of the near-inertial response in the far-field agrees 
with the excitation of an energetic first baroclinic mode. 

2. Thermocline 

Significant increases in the CW-amplitudes above the background level 
occurred in the thermocline (Fig. 4.7) at all the moorings and continued for 18-21 IP 
subsequent to storm passage. However, there are some marked differences in 
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TABLE 11 



Differences in the CW-rotating phases in the mixed 
layer between the various arrays in hurricane Frederic averaged over 5 IP. 



Arrays 


AO 

(°) 


CMA6 vs. 9 


70 


CMA2 vs. 6 


15 


CMA2 vs. 9 


63 


CMA2 vs. 3 


57 


CMA6 vs. 3 


63 


CMA3 vs. 9 


12 



thermocline current response. At CMA3, the amplitude increased to 28 cm/s as the 
storm passed through the region. After the time scale associated with the first 
baroclinic mode (tj = 25 h), the amplitude decreased by about 10 cm/s to a relative 
minima of 18 cm/s. Thus, it appears that the first mode indeed separates from the 
remainder of the modes after 25 h and there is a relative minima in the amplitudes by 
50 h (2tj). This is consistent with the linear theory of Gill (1984), who suggests that at 
2t n , the n 1 * 1 baroclinic mode will be 180 ^ out of phase with the remainder of the 
modes and the amplitudes should be near a relative minima. After 50 h (2t^ the 
amplitude began to increase slightly by about 4-5 cm/s. During the period from t2 to 
2t2 (7.5-15 IP), the CW-amplitude increases significantly to nearly 30 cm/s. Although 
the relative minima occurs at about 16.5 IP rather than the expected 15 IP, this is 
considered to be fairly good agreement between the relative extrema of the CW- 
amplitudes and the predicted time scales. There is also another peak in the amplitude 
after 19 IP that may correspond to the separation of the third baroclinic mode. The 
CW-amplitudes at CMA2 (Fig. 4.7b) decreased to pre-storm values and then increased 
again at roughly 5 IP, whereas the secondary’ maxima occurred after 7 IP at CM A3. 
One possibility is that there may be slight differences in the phase speed of the second 
baroclinic mode at CMA2 since the water is about 130 m shallower and there are some 
differences in the N profile at station 1 (Fig. 2.4). The CW-amplitude again 
approached a relative minima after 10 IP or roughly twice the time scale associated 
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' with the secondary maxima. Initially, the CW-amplitudes in the thermocline at CMA6 
and CMA9 increased to about 20 cm/s followed by a decrease of about 10 cm/s. After 
about 5-7 IP, there were similar increases in the amplitudes, which exceeded 20 cm/s at 
CMA6 and CMA9. 




Figure 4.7 CW-amplitudes (cm/s) of the near-inertial component in the thermocline at 
CMA2 (179 m, solid), CMA3 (251 m, dotted) and OTEC site (232 m, dashed) starting 
on 21 UTC 7 Sept, and ending 21 UTC 8 Oct.. 

In the thermocline (232 m) at the OTEC site, the amplitudes increased nearly 
four-fold over the background level to about 10 cm/s. However, a very large peak 
similar to that at CMA2 and CMA3 was not observed immediately following storm 
passage. The time difference in storm passage between the DeSoto Canyon and the 
OTEC site was about 1/2 IP. Taking this time differential into account, the estimated 
wavelength in the thermocline also approaches the 240 km observed by Black (1983) in 
the DeSoto Canyon region and estimated above from the horizontal phase differences 
in the mixed layer. 

C. VERTICAL STRUCTURE 

At CMA2 (Fig. 4.8a), the initial maxima were 28 and 10 cm/s in the thermocline 
and near-bottom layers as at CMA3. After 6 IP, the thermocline and near-bottom 
CW-amplitudes were equal (20 cm/s) and were roughly half the secondary maxima of 
the mixed layer amplitude. The amplitudes then gradually decreased to pre-storm 
levels after 12 IP. The thermocline and near-bottom amplitudes at CM A3 increased to 
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about 33 and 22 cm/s, respectively (Fig. 4.9b). The CW-amplitude in the near-bottom 
layers approached zero within an IP in association with an abrupt phase change. 

Within the next 4 IP, the CW-amplitude near the bottom exceeded the amplitude in 
the thermocline and then approached a maximum of 25 cm/s about 7 IP following 
storm passage. Notice that another maxima of 25 cm/s occurred in the thermocline 
after 8 IP. 

The 10 cm/s amplitudes at a depth of 950 m at the OTEC Site (Fig. 4.9c) 
demonstrates an ocean response to hurricane forcing throughout the w T ater column. 
Brooks (1983) also observed a near-inertial response close to the bottom at 700 m 
during the passage of hurricane Allen. At the OTEC site, the maximum CW- 
amplitudes in the thermocline occurred at 232 m rather than at 100 m. Because the 
CCW-amplitudes at 100 m were more energetic than the CW-amplitudes, this may be 
the influence of the shear zones on the edges of the CW-rotating Loop Current in the 
Gulf of Mexico. That is, the current meter at 100 m may have been engulfed by the 
CCW-rotating shear on the periphery of the Loop Current (Molinari and Mayer, 

1982). 

At all three moorings, the phases between the vertical levels started to converge a 
few hours prior to the closest approach of hurricane Frederic (Fig. 4.9). Generally, 
there w T as about a 90 0 phase change between the mixed layer and the thermocline at 
CMA2, CMA3, CMA6 and CMA9 as predicted by Price (1983). The near-inertial 
response below the mixed layer was nearly in-phase at all three moorings, which 
suggests a low-mode response in the vertical. Higher order modes may be present that 
cannot be detected with the limited sampling in the vertical, but the possibility of all 
the current meters simultaneously sampling at the nodes at each level is quite remote 
(J. F. Price, 1982, personal communication). The AXCP data acquired in hurricane 
Norbert to be described below provide the vertical resolution necessary to calculate the 
amount of energy in the higher order baroclinic modes. 

Since the initial phase differences were about 90° between the mixed layer and 
thermocline, the vertical wavelengths estimated from (4.2) exceeded 1000 m. Since the 
w’ater depths at the three current meter array sites were 340, 465 and 1050 m, the low- 
order modes (large vertical wavelengths) contributed significantly to the near-inertial 
variability. The 1000 m vertical wavelength is roughly five times the scale of the 
thermocline (about 200 m) and is consistent with the vertical wavelengths derived in 
Eloise simulations (Price, 1983) and the hurricane Allen observations (Brooks, 1983). 
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Figure 4.8 CW-amplitudes from a) CMA2 at 19 m (solid), 179 m (dotted), and 324 m 
(dashed) ; b) CMA3 at 21 m (solid), 251 m (dotted), 437 m (dashed), and 457 m (chain- 
dash); and c) OTEC at 100 m (solid), 232 m (dotted), 546 m (dashed), and 950 m 
(chain-dash).Time scale is the same as in Fig. 4.2. 
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Figure 4.9 CW-phascs from a) CMA2 at 19 m (square), 179 m (triangle) and 324 m 
(plus) ; b) CMA3 at 21 m (square), 251 m (triangle), 437 m (plus) and 457 m (xdash); 
and c) OTEC at 100 m (square), 232 m (triangle), 546 m (plus) and 950 m (xdash) 
every 6 h. Time scale is the same as in Fig. 4.2. 
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D. VERTICAL AND HORIZONTAL ENERGY PROPAGATION 

The vertical and horizontal group velocities at which energy propagates are 
calculated from the dispersion relationship for the baroclinic modes of near-inertial 
waves (Brooks, 1983) 

c v =-((T/m)(N 2 /<J 2 )tan 2 (<p), (4.3) 

c h = «r/K)(N 2 /<r 2 )tan 2 (tp), (4.4) 

where <7 is the observed frequency, m and k are "wavenumbers in the vertical and 
horizontal directions, and tan(cp)= k/m. The wavenumbers m and k are equal to 
2ic/X y and 2 tc/Xj 1 , respectively, where X y and represent the vertical and horizontal 
wavelengths estimated above. The mean and one standard deviation values assigned to 
the variables in (4.3) and (4,4) are : 

(7 = (7.45 ± 0.12) x 10' 5 rad/s, 

N = (6.60 ± 0.71) x 10' 3 rad/s, 
m = (6.28 ± 0.17) x 10' 3 rad/m, 
k = (2.51 ± 0.11) x 10' 5 rad/m. 

Note that in estimating the energy propagation, the Brunt-Vaisala frequency profile 
was averaged in the upper 250 m of the water column. 

The vertical group velocity (Table 12) associated with the 1000 m vertical 
wavelength is 0.15 cm/s ± 0.05 cm/s. This estimate of group velocity, which is 
associated with the first baroclinic mode, is within one standard deviation of the 0.12 
cm/s from Price (1983) and slightly larger than the estimate of 0.07 cm/s from the 
hurricane Allen observations (Brooks, 1983). The direction of the energy propagation 
is downward from the mixed layer into the thermocline for surface-intensified flow. 

For downward propagating wave groups associated with baroclinic modes, Leaman 
and Sanford (1975) showed that phase propagation is upward. After the phases began 
to separate, the vertical wavelengths were about equal to the water depth, which is 
indicative of a mode 2 response. The vertical group velocity is then 0.03 ± 0.02 cm/s, 
which is an order of magnitude less than during the hurricane forcing period. 

The horizontal group velocity of the first baroclinic mode is 0.4 m/s. The 
estimates from the Eloise simulations (Price, 1983) and Allen observations (Brooks, 
1983) were 1 and 0.3 m/s, respectively. The discrepancy between the Eloise and the 
first baroclinic mode is due to the difference in horizontal wavenumber. The rate of 
energy propagation agrees with the group velocity estimated from the Allen 
observations (Brooks, 1983). 
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TABLE 12 



Vertical (c yn ) and horizontal group velocity (c^) 
for the first two baroclinic modes estimated from the dispersion relation. 





c 

vn 


c hn 


n 


(cm/s) 


(m/s) 


1 


0.15 


0.4 


2 


0.03 


0.1 



As noted above, the rates of energy propagation depend on the horizontal and 
vertical wavenumbers (4.3, 4.4). These wavenumbers inferred from CMA6 and CMA9 
are nearly equal to those derived from the other arrays. Thus, the estimates of group 
velocity from the other two arrays are not significantly different from the above values. 

E. MODAL ANALYSIS 
1. Flat Bottom Modes 

Consider the case in which the bottom is flat, with the depth D = constant. 
The Sturm-Liouville problem is solved using a constant Brunt-Vaisala frequency (N = 
2.15 x 10'° rad/s), which was derived by vertically averaging the profile over the entire 
water column to compare to a constant N model with a sloping bottom (Lai and 
Sanford, 1986). The free mode model (Chap. Ill) is solved using a flat bottom ocean 
where the flow normal to the bottom is zero ((p n = 0 at z= -D) and a free surface 
boundary condition that pressure is continuous across the interface, 

*Ss. I k 2„ 

dz (<r 2 -f 2 ) 

Equation (3.6) along with these boundary conditions constitutes a Sturm-Liouville 
problem for which any number of modes satisfy the .relation. Since the oceanic 
response to hurricanes is dominated by the low-order modes (large vertical 
wavelength), this assumption is reasonable. Leaman and Sanford (1976) point out that 
as the vertical wavelengths (higher order modes) decrease, stratification effects become 
increasingly important. The resulting amplitudes of the horizontal velocity 
eigenfunctions for modes 0, 1 and 2 at CM A3 are shown in Fig. 4.10. The 
eigenfunctions of the horizontal velocity at all of the array sites are similar. 
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Figure 4.10 Amplitude of the horizontal velocity eigenfunctions for the flat-bottom 
mode 0 (solid), mode 1 (dashed) and mode 2 (dotted). The physical parameters are N 
= 2.15 x KpV*, c - 1.04f and D = 500 m. 

2. Bottom Slope Modes 

For the current meter measurements acquired in the DeSoto Canyon, a 
bottom slope (typically a $s = 6x10’^) is introduced into the model. Wunsch (1968, 
1969) studied the propagation of internal waves up a slope using a constant N model. 
Lai and Sanford (1986) extended the model to include rotational effects. The model 
assumes a wedge-shaped, linear bottom-slope of the form z = -a s x. The kinematical 
boundary conditions at the ocean surface and bottom in terms of the perturbation 
streamfunction (y) are: y(x,0) = 0 and y(x,-a s x) = 0 

In this case, (3.6) is no longer separable and the waves differ considerably 
from simple plane waves. Lai and Sanford (1986) showed that the propagating wave 
solutions are 

y n (x,z) = exp( ± jq(ln(yx-z)) - exp( ± jq(ln(yx + z)), 
where q = 2nn/ln (y + « s )/(y-a s ), y^ = (<x^ - f^)/(N^ - c^) is the slope of the internal 
wave characteristic, and + (-) refer to downward (upward) energy propagation. The 
perturbation streamfunction is related to the currents in the usual manner, viz. u = 
-dy/dz, w = dy/dx. The baroclinic modes are affected the most by the sloping 
bottom. 
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It should be noted that the two-dimensional, linear wedge, eigenvalue problem 
is sensitive to variations in N and x, where x > D/a . If the bottom slope is 
supercritical (a $ > y), the waves are reflected seaward and cannot propagate up the 
slope onto the shelf. A subcritical condition is appropriate at CM A3 since a $ < y. 

The resulting amplitudes associated with downward energy propagation for modes 1 
and 2 of the horizontal velocity are shown in Fig. 4.1 1. Notice the increased, near- 
bottom amplitude induced by the sloping bottom. 

Vnz 
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Figure 4.11 Amplitude of the horizontal velocity eigenfunctions for the slope mode 1 
(solid) and mode 2 (dashed). The physical parameters are the same as in Fig. 4.10, 
except that x = 85 km and a s = 6 x 10' 3 . 

3. Data Analysis 

Since the modal contributions can contain both CW and CCW-rotating 
components (Leaman, 1976), only the real part of the velocity amplitude was least- 
squares fit to the first two baroclinic, flat and sloping bottom modes after the vertical 
average was removed from the time series. The data analysis concentrates on the large 
vertical scales for two reasons : 1) only three or four observations are available in the 
column; and 2) the far-field response is dominated by low-order modes (Price, 1983 ; 
Lai and Sanford, 1986). 

The time evolution of the CW-amplitude coefficients at CMA3 indicate that 
the storm excited both vertically-averaged flow and baroclinic modes (Fig. 4.12). As 
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shown in previous studies, the baroclinic modes dominated the near-inertial response to 
hurricane Frederic and were considerably more energetic than the depth-averaged flow. 
The mode 2 current exceeded 16 cm/s and gradually decreased to 8 cm/s after about 7 
IP. The first baroclinic mode also increased significantly above pre-storm levels. The 
amplitude then increased to a secondary maximum of about 24 cm/s after 3-4 IP. 

After this secondary maxima in the mode 1 amplitude, the amplitude decreased rapidly 
as the second baroclinic mode dominated the response after 7 IP. 




Figure 4.12 Velocity coefficients at CM A3 for depth-averaged mode (solid), mode 1 
(dashed), and mode 2 (dotted). Time scale is the same as in Fig. 4.2. The biases for 
each mode have been applied to the appropriate curves to represent the amount of 
aliasing due to limited vertical sampling as described in the Appendix C. 

The depth-averaged flow reached a peak amplitude of about 7-9 cm/s during 
the spin-up phase and then decreased to a relative minimum within an IP (Fig. 4.12). 
After a few IP, the depth-averaged flow decreased to about 4 cm/s. This initial 
amplitude is comparable to the 4-8 cm/s current amplitudes induced by the free-surface 
slope, which will be addressed further in Chapter VI. It should be noted that the 
wavelength of the barotropic mode is 4700 km, whereas the basin scale in the Gulf of 
Mexico is 1500 km. Thus, it seems unlikely that the basin will admit a barotropic 
mode excited by hurricane passage. Geisler (1970) showed that the relevant parameter 
associated with the barotropic mode is the deformation radius (1000 km) rather than 
the horizontal wavelength. Although it is unclear at this point whether Frederic 
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excited a barotropic response, it certainly excited a non-zero, depth-averaged 
component. 

The modal time series are reconstructed at each depth by multiplying the 
modal coefficients by the appropriate eigenfunctions for that depth. The residual 
amplitudes are estimated by subtracting the modal series from the data (Fig. 4.13). 

The residual currents are smaller if the depth-averaged current is extracted from the 
data series prior to the least-squares fit. The residual amplitudes reach a maximum of 
7-8 cm/s after 4 IP and then decrease gradually to zero by 8 IP. 
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Figure 4.13 Residual variances for the flat bottom (dashed) and sloping -bottom 
(dotted) least squares fit with the vertical-average removed. Time scale is the same as 
in Fig. 4.2. 



The amplitudes of the baroclinic modes (Fig. 4.14) for the sloping bottom 
model were greater than in the flat-bottom case. Initially, the amplitude of the first 
baroclinic slope mode was 8 cm's larger than the flat bottom mode 1 amplitude, 
whereas the second bottom-slope mode amplitude increased by about 10 era's. The 
maximum differences between the flat and sloping bottom amplitudes occurred after 
about 2 IP and were 11 and 16 cm,'s for modes 1 and 2, respectively. The residual 
amplitudes (sec Fig. 4.13) associated with the slope mode fits generally decreased by 
about 2-3 cm's, which indicates that more of the observed variance can be explained by 
the bottom-slope modes. 



88 



10-1 



10 - 



ft- 









\\ 
\ % 



'•'*.V.?*.**j** 
-3 -1 



3 ft 

T/IP 












u 



it 



Figure 4.14 Differences in the sloping bottom and flat-bottom modes 1 (dashed), and 
mode 2 (dotted). Time scale is the same as in Fig. 4.6. 



An estimate of the amount of aliased energy was made by attempting to 
represent the currents from a 15-level baroclinic model (Adamec el ai, 1981) with only 
three currents at the same depths as at the NAVOCEANO arrays (see details in 
Appendix C). The difference in the amplitudes between the 15-level and 3-level fits is 
an indication of the bias induced by limited vertical sampling. The biases in the first 
two baroclinic modes are 6 and 4 cm/s, whereas the discrepancy between the depth- 
averages is 2 cm/s. These biases in modal amplitudes were applied by creating a time 
series with a phase of 71/2 about 5 h prior to the point of closest approach. 

The contributions of the baroclinic modes to the mixed layer amplitudes can 
be estimated from the product of the mixed layer eigenfunctions (unity in the mixed 
layer) and the modal coefficients. After accounting for the bias in the modal 
coefficients, the mode 1 and 2 contributions were about 13 and 16 cm/s, respectively. 
Since the total amplitude in the mixed layer was about 90 cm/s, the percentage of the 
amplitude accounted for by the first two modes is 33 %. Gill (1984) estimated 
amplitudes of 13 and 14 cm/s for the first two baroclinic modes lor a mixed layer depth 
of 50 m (see his Table 1). Since the initial mixed layer amplitude in that study was 100 
cm/s, the amount of baroclinic energy in the mixed layer was 27 % of the total. 
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During the period of the enhanced mode 1 amplitude (3 IP), the total contribution of 
the baroclinic modes to the mixed layer was about 37 cm/s or roughly 41 %. Thus, the 
amplitudes associated with the first two baroclinic modes estimated from Frederic data 
are within 10 -15 % of theoretical results derived from linear theory. 

The modal coefficients are then used to estimate the amount of variance 
associated with the summation of the depth-averaged flow and the first two baroclinic 
modes. The evolutions of the explained variances averaged over an IP are given in 
Table 13 for the flat- and sloping-bottom modes at CM A3. The flat-bottom mode 1 
contributed 38% while the bottom-slope mode 1 accounted for 41 % of the variance 
over the first IP. The contributions from the first mode increased to a maximum over 
the following 3-4 IP for both the flat and sloping bottom cases. There was also an 
increase in the amount of explained variance by mode 2 (30-40 %) after 7 IP. Similar 
contributions from mode 2 are found in the sloping-bottom model with a further 
increase in the explained variances by 6 %. Hence, the increases in the thermocline 
amplitude 6-8 IP after storm passage appear to be associated with the second 
baroclinic mode. 



TABLE 13 

Explained current variances at CMA3 for the depth-averaged and first two 
baroclinic modes in the flat-(F) and sloping-(S) bottom models. 







F 






S 


F 


S 


IP 


0 


1 


£(1 + 2) 


1 


£(1 + 2) 


L( 0 + 


1 + 2) 


0 


0.33 


0.50 


0.86 


0.49 


0.87 


0.80 


0.83 


1 


0.34 


0.38 


0.72 


0.41 


0.80 


0.65 


0.69 


2 


0.18 


0.40 


0.79 


0.36 


0.S4 


0.72 


0.78 


3 


0.27 


0.60 


0.88 


0.53 


0.89 


0.80 


0.85 


4 


0.19 


0.68 


0.83 


0.65 


0.82 


0.81 


0.82 


5 


0.13 


0.62 


0.69 


0.73 


0.66 


0.71 


0.68 


6 


0.06 


0.54 


0.55 


0.69 


0.59 


0.57 


0.59 


7 


0.19 


0.24 


0.42 


0.46 


0.54 


0.31 


0.45 


8 


0.20 


0.11 


0.56 


0.11 


0.60 


0.43 


0.55 


9 


0.24 


0.20 


0.50 


0.15 


0.64 


0.47 


0.60 


10 


0.23 


0.17 


0.45 


0.25 


0.51 


0.57 


0.62 
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A sum of the depth-averaged flow and first two baroclimc modes at CM A3 
explained 62 % and 68 % of the variance over the first 10 IP for the flat- and sloping- 
bottom models, respectively. The depth-averaged current contributed about 21 % to 
the observed near-inertial wave variance over the first 10 IP. The first baroclinic mode 
accounted for 40 % of the current variations, whereas the second baroclinic mode 
contributed about 25 % in the flat-bottom case. By contrast, the first baroclinic 

bottom-slope mode accounted for about 45 % of the near-inertial variance. Thus, the 

✓ 

increases in the explained variances by the first two bottom-slope modes were about 5 
% and 6 % , respectively. 

At CMA2, the depth-averaged term accounted for 25% of the variability of 
the near-inertial response, and a sum of the modes explained 52 % of the variance. 

The second baroclinic mode was more energetic than mode 1 by about a factor of 3. A 
sum of the modes at the OTEC site contributed about 62 % to the near-inertial 
variance. However, some caution has to be applied to these results in light of the 
intermittency of the forced near-inertial response at 100 m. 

At CMA6 and CMA9, the depth-averaged components ranged from 5-9 cm/s 
and accounted for about 20 % of the observed near-inertial variance. The amplitudes 
of the first two baroclinic modes were within 3-5 cm/s of those estimated from the 
modal analyses at CMA2 and CMA3. Furthermore, the modal response was 
dominated by the first baroclinic mode. 

In comparison, the modal decomposition of vertical current profiles into the 
first three baroclinic modes during the passage of hurricanes Dawn and Carrie 
explained 55 and 73 % of the variance for the flat and bottom-slope models, 
respectively (Lai and Sanford, 1986). Their results clearly showed a fairly significant 
change in modal contributions by including a sloping bottom. The velocity profiles 
collected during STREX show that higher order modes are excited by the passage of a 
storm, but most of the energy seems to be concentrated in a mode 2 response 
(D'Asaro, 1984). The low-mode response determined from the Frederic observations is 
consistent with these earlier studies except that the depth-averaged mode also 
contributes to the variability. 

F. SUMMARY 

Analyses of the moored current meter measurements indicate that the forced 
near-inertial waves have frequencies that are blue-shifted above the local inertial 
period, large CW-rotating amplitudes and wavelengths of the first baroclinic mode. 
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This estimate of wavelength (250 km) is similar to the findings of Black (1983) using 
the AXBT measurements made in the vicinity of the current meter arrays. 

The time-evolution of the first two baroclinic modes is similar to the predictions 
from linear theory (Gill, 1984). The secondary increases in the CW-rotating amplitudes 
in the mixed layer and in the thermocline are correlated to the time scales associated 
with the first two baroclinic modes. The vertical structure of the near-inertial response 
can be described by a sum of the depth-averaged and first two baroclinic modes. The 
addition of the bottom slope accounts for 5-8 % more of the near-inertial variability. 

The critical issue raised by these calculations is the role of the depth-averaged or 
barotropic component in the near-inertial frequency band. At all of the arrays 
(including OTEC), there was a significant increase in the current amplitudes during the 
passage of Frederic. Because of limited vertical sampling, some higher mode energy 
may have been leaked into the barotropic mode which has amplitudes of 7-10 cm/s. 
This issue is addressed in Chapter VI using numerical simulations from a free surface 
model. 
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V. VERTICAL STRUCTURE 



The primary motivation for this chapter is to determine the number of baroclinic 
modes that are required to describe the near-inertial ocean current response in the 
near-field. Since the current meter moorings only provided data at selected depths, the 
least-squares fit may have aliased some of the higher mode energy into the low 
baroclinic modes. Thus, the near-inertial response in the ocean current profiles 
acquired during hurricane Norbert are decomposed into the first three or four 
baroclinic modes using the forced, linear theory in Chapter III. 

A. SYNOPTIC SCALE PICTURE 
1. Mixed Layer 

Previous studies of the oceanic response to hurricane passage using AXBT 
observations have documented a crescent-shaped pattern of sea-surface temperature 
patterns (Black, 1983). The maximum temperature decreases of 2-3 °C in the mixed 
layer occur on the right side of the hurricane at about 1-2 R^^. behind the storm 
center. Numerical studies also demonstrate this rightward bias of upper ocean 
processes (Chang and Anthes, 1978 ; Price, 1981). 

In hurricane Norbert (Fig. 5.1a), the maximum cooling region is located to the 
rear of the storm at roughly 2R max . The minimum mixed layer temperature is 25.7 
°C, which represents a decrease of 2.6 °C relative to the observed temperature ahead of 
the storm. Over the remainder of the domain, the mixed layer temperatures exceeded 
27 °C. 

The distribution of mixed layer depths in Norbert correspond to the pattern of 
mixed layer temperature decreases (Fig. 5.2). The mixed layer depth decreased only to 
18 m at roughly 90 km (2-4 R^^) in the wake of the storm where peak upwelling 
occurs. Strong vertical advective processes are associated with the upwelling of cool 
water, which tends to form a shallower mixed layer and indirectly enhances 
entrainment mixing and surface cooling. On the right side of the track, the mixed layer 
deepened to depths of 40 m or more. 

Maximum mixed layer currents in the wake of Norbert exceeded 1 m/s at 
roughly 2-4 R max (Fig. 5.2b) in the same region as the minimum in mixed layer 
temperature and depth. Sanford et al. (1987) also observed similar maximum velocities 
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primary and secondary radii of maximum winds as depicted by the quasi-circular 
dashed lines and the direction of the storm movement is indicated by the double arrow. 



in the three-layer model fits to the observations. The general pattern of mixed layer 
currents agrees with the divergent upwelling and convergent downwelling regimes 
predicted from numerical studies (Price, 1981) and shown in analytical studies (Geisler, 
1970). Undoubtedly, larger current speeds existed at the surface, but these could not 
be measured by the AXCP. 

2. Thermocline 

The thermocline current (80 m) pattern ii N 'rbert was anticyclonic and 
exceeded 40 cm/s (Figure 5.3). In the region of maximum mixed layer currents and 
strong upwelling, the difference in the directions of the currents between the mixed 
layer and thermocline is only 40 -60 °. Over the remainder of the domain, the direction 
of the thermocline currents is nearly opposite to the mixed layer current, which 
suggests a low baroclinic mode response (discussed below). The considerable vertical 
current shears in the upper ocean are presumably associated with near-inertial 
processes. Thus, a fairly large area of near-critical Richardson numbers (shown below) 
( < 1) is located in the regions of largest mixed layer temperature decreases. These 
calculations of Richardson numbers are the first confirmation of this physical process 
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Figure 5.2 Mixed layer depths (m) (left) and mixed layer currents (cm/s) in Norbert 
with the shaded regions depicting the regions of maximum mixed layer depths (> 40 
m) and currents ( > 120 cm's). 



that has been generally assumed in numerical treatments of the oceanic response to 
hurricanes (Price, 1981). However, surface-generated mixing is also a possibile 
explanation because of high surface winds and very shallow mixed layers. 

B. VERTICAL STRUCTURE 

1. Vertical Wavenumber Spectra 

A Tukey data window (Otnes and Enochson, 1978) is applied to the current 
profiles prior to the removal of the mean in the calculation of the energy spectra. The 
data are transformed and spectrally-averaged over bandwidths to help decrease leakage 
of energy to adjacent frequency bands (Otnes and Enochson, 1978). Energy above the 
Nyquist wavenumber (l/2Az where Az is the 3-m sampling interval) is eliminated to 
minimize aliasing of the spectra. 

The vertical wavenumber (m) spectra of the kinetic energy (KE) are compared 
to the background internal wave spectra of Garrett and Munk (1975). The form of the 
spectra used in this intercomparison is similar to that of Leaman (1976) 

KE(m) = (N 0 2 /M 3 )(3E 0 /2j,7t)A(X v ), (5.1) 



95 




Figure 5.3 Currents (cm/s) at 80 m in hurricane Norbert •with the shaded area depicting 
the region of maximum currents ( > 40 cm/ s). 

where A(X y ) = (t-l)(H-X v )‘\ = m/(j*7tM o ), M o = 0.122 cycles/km, j*= 6, E Q = 

In x 10'^ and and t= 2.5 (Garrett and Munk, 1975). The nondimensional energy 
density (E q ) seems to be a universal value within a factor of 2-3 in the ocean for the 
background internal wave spectra (Munk, 1981). 

In the right-front quadrant of hurricane Norbert (AXCP 02 in Fig. 2.6), the 
observed vertical wavenumber spectra are not significantly diiTerent from the GM75 
spectra in the wave band from 50 m to 350 m (Fig. 5.4). At the large and small 
wavenumbers, there are only small differences between the observed and GM75 
spectra. However, the observed KE spectra in the left-rear quadrant (AXCP 20), 
exceeds the GM75 values by nearly an order of magnitude over the entire spectrum. 

The most energetic vertical wavenumbers are associated with vertical wavelengths 
greater than 50 m. At the large vertical wavelengths (> 500 m), the KE values are 
more than an decade greater in the wake than in front of the storm because of the 
increased internal wave activity induced by the moving hurricane (Geisler, 1970). 

The energy spectra are decomposed into clockwise (CW) and counter- 
clockwise (CCW) rotating components (Mooers, 1973). In the front of hurricane 
Norbert (AXCP 02), the CW- and CCW-rotating components are equal in magnitude 
(Fig. 5.5a). In the wake of Norbert (Fig. 5.5b), the CW-rotating energy exceeds the 
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Figure 5.4 WKBJ-sti etched vertical wavenumber (cpm) spectra of KE (cm/s) 2 /cpm 
(solid) in front of hurricane Norbert AXCP 02 (loft) and in the wake of hurricane 
Norbcrt AXCP 20 (right). I he Garrett and Munk (1975) spectrum is represented by 
the dashed line. The 95 % confidence limits are also shown. 



CCW-rotating by an order of magnitude in the range of wavelength range of 50-200 m. 
At larger ( > 500 m) and smaller ( < 50 m) wavelengths, there is no preference of 
rotation as the CW and CCW rotating energies are equal. 

2. Mean Current Profile 

The spatial-mean currents are estimated from the Norbert AXCP's to separate 
the effects of the hurricane from the background current variability or pre-storm flow 
conditions. Prior to estimating the spatial mean, the vertical average is removed from 
each profile (Leaman and Sanford, 1975). Bootstrap methods are used to calculate the 
mean profiles of currents and Brunt-Vaisala frequency profiles in the uppper levels 
until less than nine profiles are available (Effrom and Gong, 1983). To estimate a 
deeper mean profiles, the mean profiles are calculated using arithmetic means until less 
than three profiles are available for averaging. This averaging process is concluded at 
660 m. The uncertainities in estimating the mean are calculated by assuming a normal 
distribution of the mean currents. 

The two-layer flow regime in the mean currents in hurricane Norbert (Fig. 5.6) 
is fairly consistent with the expected climatology north of the equatorial currents in the 
Pacific. Maximum mixed layer currents are about 16 cm/s towards the north-northeast 
with current reversals at about 100 m in the v-component and at 60 m for the u- 
component. Both mean current components approach zero as depth increases. 

3. Near-Inertial Current Profiles 

To isolate the near-inertial (and higher frequency) currents, the spatial-mean 
current profiles are subtracted from each of the profiles (Fig. 5.7). In front of the 
storm at r = 2R max (AXCP 31), the wind stress produces a u-component of 25 cm/s 
towards the west. There is considerable vertical structure in the v-component with 
vertical wavelengths of about 200 m superposed on the larger vertical wavelength 
features. In the center of the storm (AXCP 15), the wind is forcing a very large u- 
component towards the west. Since this probe descended slowly, the profile only 
extends to about 240 m. The vertical structure of the u-component is fairly simple 
compared to that associated with the v-component. Since AXCP 04 was deployed at 
roughly r = 6 R max in front of the storm, the wind stress and the ocean velocities are 
considerably weaker than those observed closer to the center of the Norbert. There are 
fairly strong currents in the upper ocean ( < 100 m). The maximum upper ocean u and 
v components of currents are 20 and 40 cm/s directed towards the northwest. Below 
200 m, the current components are nearly equal in amplitude (10-20 cm/s) with shorter 
vertical scale oscillations (100-200 m) in the v-component. 
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Figure 5.5 WKBJ-streclied vertical wavenumber (cpm) spectra of the CW- (solid) and 
CCW-rotating (dashed) KE (cm 2 /cpm) spectra (cpm) for hurricanes Norbcrt AXCP 02 
(left) and in the wake of hurricane Norbcrt at AXCP 20 (right). The 95 % confidence 
limits are also shown as in 5.4. 




Figure 5.6 Spatially-averaged u-component and v-component (cm/s) for hurricane 
Norbert. The 95 % confidence limits are based on a Gaussian distribution and 
bootstrap methods. 



One of the features predicted by numerical models is that maximum currents 
are in the region of r = 2R max to the right of the storm track (Chang and Anthcs, 
1978). In the right-rear quadrant of Norbert at r = 3R max (AXCP 14), the maximum 
horizontal current velocities are about 60 cm/s and suggest a simple vertical structure. 
In the wake of the storm in the left-rear quadrant (AXCP 20), the v-component 
increases to 40 cm/s towards the south in the upper ocean, whereas the u-component 
exceeds 75 cm/s and is flowing away from the track. This velocity pattern is indicative 
of divergent motion and upwelling processes. The velocities exceed 20 cm/s below 400 
m. which may indicate vertical energy propagation from the mixed layer. 

4. Richardson Number Profiles 

In numerical models of the oceanic response to hurricanes, mixing is 
sometimes parameterized in terms of a Richardson number (Price, 1981). If the 
Richardson number falls below a certain critical value, momentum and heat are mixed 
between the two layers. Local gradient Richardson number profiles in Norbert can be 
calculated from the near-inertial profiles described in Fig. 5.7 and the estimated Brunt- 
VfiisfilS frequency profiles. The gradient Richardson number is given by 
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Figure 5.7 The near-inertial u-component (top) and v-component (lower) profiles 
(cm/s) from hurricane Norbert at selected locations as indicated in Fig. 2.6. 
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where u z and v z represent the vertical shear of the near-inertial profiles in Fig. 5.7. 

Since large-scale, near-inertial waves are addressed here, a Az of 9 m is chosen. 

D'Asaro (1985) argues that this interval represents a compromise between the inherent 
error in the measurement of velocity shear and temperature and an optimal scale. 

Profiles of Richardson number (Fig. 5.8) tend to be erratic because of the 
possiblity of intermittent shear instabilities (Desaubies and Smith, 1982). For example, 
layers in which the gradient Richardson falls below a critical value of 1/4 suggests that 
conditions are favorable for shear instabilities. At AXCP 31 (in front of the storm), 
shear instability is enhanced between 100 and 200 m, which corresponds to a large 
change in the vertical temperature gradient from 0.12 °C/m to 0.07 °C/m. At AXCP 
15 (center of the storm), the Richardson numbers fall below the critical value in the 
layer between 20 to 60 m and 120 m. Since the mixed layer depth at AXCP 15 is 
greater than 40 m (see Fig. 5.2), shear instability may be one of the mechanisms 
responsible for mixed layer deepening. At AXCP 04 (right-front quadrant), the 
Richardson number falls below the critical value intermittently between 200-300 m 
which suggests some small scale mixing events. The profiles of Richardson numbers 
are virtually identical in the upper 80-100 m at AXCP 15 and AXCP 14 . In the area 
of strong vertical advection (upwelling) at AXCP 20, the shear instability is apparently 
occurring throughout the thermocline. This region is still being strongly forced because 
the winds are high in a region of shallow mixed layers. 

The near-inertial velocity profiles acquired during the passage Norbert indicate 
the presence of large vertical wavelength oscillations. Since these profiles were 
de A \C" ed within the direct forcing regime or near-field, it is clear that vertical structure 
is not just a sum of the free baroclinic modes as in the far-field (Gill, 1984). The wind 
stress is a very' potent source of near-inertial oscillations in the upper layers (D'Asaro, 
1984). Thus, the current profiles should be expanded in terms of the applied 
atmospheric forcing, especially for comparisons in the near-field. 

C. MODAL ANALYSES 

Gill (1984) noted that if the large scale limit (S> 1) is satisfied (see Table 4), the 
ocean response can be described by the first few free modes. As shown above, 55-65 
% of the variance in the current meter observations in the far-field after the passage of 
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Figure 5.8 Richardson number profiles from hurricane Norbert based on the vertical 
shear of the profiles in Fig. 5.7 and the Brunt- Vaisala frequency profile at each location 
(as indicated on Fig. 2.6). 



hurricane Frederic could be explained by a summation of the depth-averaged and first 

two baroclinic modes. However, the observations in Norbert are clearly in the near- 

field. Furthermore, mixed layer currents (N(z) = 0) are set to zero when scaled 

according to WKBJ theory, which eliminates the most energetic current components. 

Thus, the near-inertial current velocities from the AXCP's should be compared to the 

forced baroclinic modes. The residual current (u r ) at each depth that cannot be 

explained by the forced model is given by: 

u r (x\y\z)= u 0 (x\y',z)- Yu n (x',y')<p nz (z) , (5.3) 

where u o is the observed u-component rotated 40 ° clockwise into storm coordinates 

(x ,v') and u is the horizontal structure coeflicient as a function of mode number n. 
n 

The residual v-component has a similar expression. The vertical average is removed 
from the observations prior to the fitting procedures. 

For each mode, the horizontal structure functions of velocity have three terms in 
(3.13) : the wind stress curl, divergence and Ekman-type terms. In the superposition 
process, the wind stress is divided by these equivalent forcing depths. Since the 



equivalent depths (J^ nz 2 U)dz) for the first five baroclinic modes are 98, 334, 336. 421 



and 592 m (see Table 6), the wind stress and the simulated horizontal velocity structure 
functions decrease as mode number increases. 
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1. Forced Modes 

Any combination of the three forcing terms for each velocity can be used to 
construct a model profile at (x',y') using (3.9). Since the model is linear, the same 
terms for each mode must be used in the summation over n modes. Selected velocity 
profiles in the wake are compared in Fig. 5.9 to the analytical model profiles using 
three or four modes. These horizontal structure functions are based on physical 
processes calculated in (3.19), rather than as a least-square fit of orthogonal 
polynomials. 

2. Vertical Profiles 

On the right side of the track (AXCP 13 , see Fig. 2.6 for locations), all three 
forcing terms are used for both velocity components in a summation that included the 
first four baroclinic modes (Fig. 5.9a). The model profile can account for over 70 % 
of the observed variance using the first four near-inertial modes! In the near-surface 
layer (20-40 m), the observed velocities exceed the model velocities by only 15 cm/s. In 
the 50-70 m layer, there is a difference of 15-18 cm/s between the observed and model 
v-component, which suggests that higher order near-inertial baroclinic modes may be 
present, or that all of the surface wave energy has not been removed. The enhanced 
vertical shear between 80-100 m is associated with near-inertial processes. Within 10% 
uncertainity limits of the model storm parameters, the residuals in the upper layers are 
reduced below 15 cm/s. The depth-average residual current is 4 and 3 cm/s for the u 
and the v-components, which is indicative of a reasonably good fit. On the right side 
at AXCP 14 (Fig. 5.9b), a four baroclinic mode model also agrees a little better with 
the observed velocities in the enhanced vertical shear layer (50-70 m). 

On the left side of the track at r = 2R max (AXCP 18), there is considerable 
agreement between a Tnee mode model and the observed profiles (Fig. 5.9c). The v- 
component driven by the wind stress curl changes direction between the upwelling and 
downwelling peaks in the wake. The model u-component in the upper 100 m 
underestimates the observed profile by 18-20 cm/s although the depth-averaged residual 
is about 5-6 cm's. Compared to the right side, the agreement between the observed 
and simulated profiles at this site is poor. Since hurricane Norbert changed both speed 
and direction (curving towards the west) in the vicinity of the AXCP deployments (Fig. 
2.6), the oceanic response on the left side of the storm will be considerably more 
complex than on the right side. 
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Figure 5.9 Observed (solid) and the linear model profiles (dashed) in cm/s for the u- 
component (upper) and the v-componcnt (lower) for selected AXCF's deployed in the 
wake of hurricane Norbert. The model profiles are based on a summation of three or 
four baroclinic modes as discussed in the text. 



In the left rear-quadrant at r = 3R max (AXCP 20), the model profile with 
three modes overestimates the observed u-component in the mixed layer by 5-10 cm/s 
and underestimates the v-component by about 20 cm/s (Fig. 5.9d). Below the mixed 
layer, the residual velocities are within 5-10 cm/s and are represented reasonably well 
■with the near-inertial model. However, there is a large discrepancy at 180 m between 
the model and observed profiles. Adding more vertical modes to the summation does 
not significantly reduce the residual profiles. In this region (see Figs. 5. 1-5.3), there is 
a minimum in mixed layer temperatures and depth (18-24 m) as upwelled water from 
the thermocline is entrained into the mixed layer. The vertical velocity of the first 
mode (see Fig. 3.3) is positive, which is associated with an upwelling regime and the 
divergence in the upper ocean (relative maxima of u-component). 

On the left side at r = 3-4 R max (AXCP 21), considerable vertical structure 
exists in the observed u-component between 80 and 120 m that is not well resolved by 
the linear model using four baroclinic modes (Fig. 5.9e). The model u-component is 
within 8-10 cm/s of the observed profile. However, the mixed layer v-component is 
underestimated by over 30 cm/s. Slightly farther back in the wake (r = 5-6 R max )> the 
model v-components are comparable to the observed v-component in the mixed layer 
(not shown). One possible explanation for the large discrepancy in this region could be 
associated with the wind shifts associated with the curved path of the hurricane (see 
Fig. 2.6), which should affect the profiles more on the left side than on the right side. 
Nevertheless, the averaged residuals are 4-5 cm/s for both profiles. The shorter vertical 
wavelengths (200 m) in the observed profiles cannot be resolved with only four 
baroclinic modes. However, these smaller scale features are not nearly as important in 
terms of energetics as the large vertical wavelength modes. 

A three mode model follows the observed velocity profile on the left side of 
the storm at AXCP 03, except for the 200 m vertical wavelength oscillations (Fig. 5.9f). 
The residual u-component is less than 10 cm/s throughout the profile with an average 
residual of 4 cm/s using the three modes. Again, the very large v-component near the 
surface is underestimated by about 40 cm/s and there are fairly large residuals 
throughout the upper 100 m. 

3. Kinetic Energy 

The observed kinetic energy at selected AXCP sites is compared in Table 14 to 
the simulated energies for the model with three or four modes as discussed above. The 
most striking feature is that mode 1 contributes more to explaining the observed energy 
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than any other mode. On the right side of the storm, the model energy explains 
between 67-78 % of the observed energy at AXCP's 13 and 14 when the first four 
baroclinic modes and all of the forcing terms are used in the model. Thus, the stress 
divergence and Ekman-type terms contribute to the ocean current variability, but they 
are dominated by the wind stress curl terms in (3.19). A least-squares fit of the 
observed velocity profiles at AXCP 13 and 14 to the first three modes accounts for 80 
and 55 % of the observed KE. On the left side of the storm, there is considerable 
spatial variability in observed KE, and the model accounts for 34-75 % of this energy. 
The minimum in the observed KE at AXCP 18 suggests that it was deployed in a 
region in which the v-component changes sign. The analytical model profiles with 
three baroclinic modes only explains 34 % of the observed kinetic energy as compared 
to 20 % in a least-squares fit. A least-squares fit of the first ten modes only explains 
55 % of the observed kinetic energy. At AXCP 20 and 21, the analytical model 
explains over 60 % of the observed kinetic energy using three baroclinic modes. Thus, 
it is clear that the oceanic response changes dramatically over relatively short spatial 
distances (1-2 R max )- At AXCP 03, the model kinetic energy only explains half of the 
observed energy using either a three or four mode model. At this location, the least- 
squares fit explains about 65 % of the energy because the mode 1 coefficients are 
about 10-12 cm/s more than the simulations from the analytical model. Some of the 
discrepancies between the observed and analytical profiles are due to inherent spatial 
variability of the hurricane. For example, Norbert was changing both speed and 
direction just south of the vicinity where the AXCP's were deployed. For a hurricane 
curving towards the west (Fig. 2.6), the oceanic response on the left side will be 
considerably more complex. Furthermore, the actual wind stress fields were not 
available mainly due to limited observations in the left-front quadrants. A more 
quantitative comparison between the observed and model profiles may require that the 
Rankine vortex be replaced with the actual spatial distribution of the wind stress. 

4. Model Sensitivity 

To assess the impact of uncertainties in the observed storm parameters on the 
model simulations, the model forcing functions (t , R max > U h , and inflow angles) are 
changed by ± 10 % for the first three baroclinic modes (Table 15). A composite of the 
changes in the simulated velocities for a 10% change in the storm parameters is 
combined for each baroclinic mode and then summed over the first three baroclinic 
modes. These changes in model parameters are combined into a composite value for 
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simulations 


TABLE 14 

Percentage of the observed KE explained by the 
based on a summation of three or four baroclinic modes at 
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each mode. As expected, the differences in the simulated velocities are largest closest 
to the storm (3-7 cm/s) and are dominated by mode 1. Surprisingly, the uncertainties 
in the velocity for mode 3 exceed that in mode 2. Since the wavelength of mode 3 in 
the unforced case is about 78 km (compared to 130 km for the second baroclinic 
mode), the mode 3 wavelength is closer to the scale of the wind stress curl, which may 
account for the increased excitation of the modal amplitude. 



TABLE 15 



Change in the horizontal structure functions (cm/s) 
for the first three baroclinic modes based a 10 % uncertainty in storm 

parameters. 
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D. SUMMARY 



The near-inertial response is isolated from the current profiles by subtracting the 
spatially-averaged mean currents. The frequency of the waves was determined using a 
mixed layer Burger number as in hurricane Frederic. The frequency of the near-inertial 
waves (<r= 1.180 is within the internal wave band of 0.9f-1.2f (Mooers, 1975; Kunze, 
1985). 

On the right side of the storm track, the forced, linear model explains over 70% 
of the observed oceanic current variability using th6 first four baroclinic near-inertial 
modes. The enhanced vertical shear between 50 -100 m appears to be associated with 
near-inertial wave processes. Some of the large differences between the observed and 
model velocities are due to uncertainities in the wind stress forcing that is used to 
represent the hurricane. Since the hurricane is modelled using a Rankine vortex, all the 
complexities and spatial variations in the wind stress are not reproduced. All of the 
terms (wind stress curl, stress divergence and Ekman-type) are used to explain the 
currents measured by the AXCP's 13 and 14. 

On the left side of the storm track, the observed profiles are generally best fit 
with a model based on only three baroclinic modes. The mixed layer v-components on 
the left side are underestimated by 15-30 cm/s, whereas the u-components are within 10 
cm/s of the observed. Since the path of hurricane Norbert curves to the left in the 
vicinity of the AXCP deployments, the ocean response of the left side is considerably 
more complicated than assumed in the model physics. The observed current profiles at 
AXCP 18 are particularly complicated because this position is in a strong upwelling 
and mixing regime. In a qualitative sense, the simulated v-component agrees with the 
observed profile. The residuals at both AXCP 18 and 20 are 4-6 cm/s and the 
enhanced shear layer is modelled fairly well at AXCP 20. Farther away from the storm 
(r = 3-4R max ), the model u-component nearly equals the observed velocity, even 
though the observed v-components are still 20-40 cm/s larger than predicted from the 
model. Below 100 m, the observed velocities are modelled fairly well with averaged 
residuals of 3-4 cm/s. Consequently, it appears more complicated dynamics are 
involved on the left side of Norbert, especially in the region of maximum shallowing of 
the mixed layer. As noted in Chapter II, the Rossby number for the Norbert storm 
was 0(0.5) which implies that the response may be moderately nonlinear. In 
particular, mixing and nonlinear effects that have been neglected in this simple 
treatment may be important near the center. 
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VI. NUMERICAL MODELING 



Chang and Anthes (1978), Adamec et al. (1981), Hopkins (1982), Price (1983) 
and Greatbatch (1983) numerically modelled the baroclinic response to a translating 
hurricane by imposing a rigid lid at the sea-surface. A deep ocean was assumed in all 
cases except in the study by Hopkins (1982), who considered an ocean depth of 964 m. 
Hopkins simulated the dynamic response to a model hurricane moving at the same 
translation speed as Frederic. The simulated mixed layer currents agreed well with the 
observations. However, the computed magnitudes of the currents in the thermocline 
were too small. It is believed that one source of the discrepancies between the 
simulations and the observations is due to the neglect of the free surface and the 
barotropic mode, which is shown to be important in the Frederic moored current meter 
observations. 

The primary motivation for using a numerical model here is to understand the 
role of the depth-averaged and baroclinic components in the ocean current response to 
hurricane passage. In addition, the numerical simulations provide a more complete 
(four-dimensional) picture of the hurricane-induced response than can only be inferred 
from sparse measurements. Because of the limited vertical sampling by the moored 
current meter arrays in Frederic, some of the higher mode energy may have been 
aliased into the estimates of the depth-averaged and first two baroclinic modes. The 
near-inertial components are isolated by demodulating the ocean current simulations 
and are fitted to the first few vertical modes. Many other interesting comparisons 
between these model simulations ard the observations are left for future studies. 

A. MODEL DESCRIPTION 

A three-dimensional nonlinear, primitive equation model with a free surface and 
flat bottom has been developed by Chang (1985) to simulate the ocean's baroclinic and 
barotropic response to hurricane passage. The ocean is also assumed to be in 
hydrostatic balance, incompressible, and on an f-planc at 29 °N. The set of governing 
equations can be found in Chang (1985). 

1. Model Grid and Storm 

A staggered Arakawa-C grid (Haltincr and Williams, 1980) is used in the 
numerical model with a Ax of 20 km. The C grid has been shown to produce realistic 
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phase and group velocities in the geostrophic adjustment process (Schoenstadt, 1977). 
The model domain is 4000 x 1000 km with cyclic boundary conditions on the east and 
west sides and Neumann conditions on the north and south boundaries. A vertical 
grid with 17 levels is stretched for increased resolution in the upper ocean (Az = 10 m) 
and a Az = 80 m at depths below the thermocline. The total depth of the ocean is 605 
m and the drag in the bottom boundary layer is treated using Rayleigh friction terms. 
The value of the friction is 2.3 x 10'^, which corresponds to an e-folding scale of five 
days. 



The model storm is represented by a Rankine vortex using the parameters 
observed in Frederic (Pow'ell, 1982) which moves westward at 6.5 m/s. Since the inflow 
angle is 20 °, there can be a significant wind stress divergence (Mayer et al . , 1981). 

2. Mode Splitting 

Since the model includes a free surface, a barotropic mode will be excited by 
the model hurricane. Madala and Piasek (1977) treated the primitive equations in a 
free-surface model using a semi-emplicit formulation. By contrast, Chang (1985) 
introduced a mode-splitting technique by subtracting the weight of the fluid 
(hydrostatic approximation) from the total pressure to form a relative pressure p r (z). 
This term is decomposed into expressions for the mean relative and perturbation 
pressures: 



Pr = 



1 

2 



gn + g £z + p eq 



Pr = Pr * Pr 



( 6 . 1 ) 



where = 1/Dj( ) dz represents a vertical average, £ is the density deviation (p-p o ), 
i] s is the perturbation height of the ocean surface from a mean depth D and p s is the 
ocean density near the surface. 

The barotropic mode is treated in the momentum balance by separating the 
equations into a mean and perturbation quantities. The vertical mean pressure 
gradient excites a barotropic mode. Accordingly, the spatial gradients of the free 
surface height induce a "vertically-averaged" or a barotropic mode into the simulations. 
The perturbations of the relative pressure excite only the baroclinic part of the ocean 
response (Chang, 1985). Since the free-surface height enters into the pressure terms, a 
prognostic equation for i] s has to be derived by assuming that w = 0 at z= -D and 
vertically integrating the conservation of mass equation 
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( 6 . 2 ) 
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The inclusion of the free surface in the ocean model requires a special time 
integration treatment because the speed of the barotropic mode is 100 times faster than 
the first baroclinic mode. To maintain numerical stability, the Courant-Frederichs- 
Lewy condition is 
CAt 

f < V2, (6.3) 

Ax 

where c fi is the phase speed of the waves to be resolved in the model. For a barotropic 
mode in 600 m of water, c q = 75 m/s, whereas the phase speed of the first baroclinic 
mode is typically 1-3 m/s. For a grid spacing (Ax) of 20 km, the maximum time steps 
would be 267 and 6667 sec for the barotropic and first baroclinic modes, respectively. 
Hence, for every baroclinic mode time step, 25 time steps are required to integrate the 
barotropic mode so that the modes can interact. 

The model is integrated in time using a leapfrog scheme and centered 
differencing in space with second order accuracy (Chang, 1984). As described above, 
there are Wo time steps: a 40 s time step for the barotropic integrations; and 1200 s for 
the baroclinic mode integration. A horizontal diffusion coefficient of 10 5 cm 2 /s is 
applied to smooth the solutions. The total time integration in the following 
simulations is 132 h or roughly 5 IP. 

3. Mixing Effects 

In upper ocean response to hurricane studies, various types of mixing 
parameterizations have been proposed. Mixing effects were shown to be important in 
the oceanic thermal response to the passage of a hurricane (Elsberry et al., 1976). 
Whereas O'Brien (1967) used K-theory (or first order closure) to close the system of 
equations, Elsberry et al. (1976) parameterized the net generation of turbulent kinetic 
energy in terms of the third power of the friction velocity (u* 3 ) times an exponential 
decay factor that was proportional to the mixed layer depth. These studies revealed 
that advective processes dominate near the center of the storm, but mixing effects are 
more pronounced outside R max * Chang and Anthes (1978) also parameterized mixing 
in terms of u* 3 . Adamec et al. (1981) embedded the mixed layer model of Garwood 
(1977) into the ocean general circulation model (OGCM) of Haney (1974). This study 
also demonstrated that mixing effects dominated advective effects outside R max - In 
contrast, Price (1981) modeled the sea-surface temperature (SST) response to a moving 
hurricane using a parameterization in terms of velocity shear (6v) at the base of the 
mixed layer based on a bulk Richardson number approach (Pollard et al ., 1973). The 
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upward heat flux (ocean losing heat to the atmosphere) was only 15% of the 
entrainment heat flux (w'T) h at the base of the mixed layer. Greatbatch (1984) 
showed that the surface sensible and latent heat fluxes were only 22% of the 



In the present version of the Chang (1985) model, vertical mixing is 
parameterized using an eddy diffusivity K z based on mixing length (l m ) and the local 
Richardson number (Ri): 



where the gradient Richardson number is defined in (5.2). The mixing length is 
assumed to be the depth of the mixed layer, which thus represents a rather large eddy 
scale. The vertical fluxes at the top and bottom of the ocean are assumed to be zero in 
this model. 

4. Initial Stratification 

The model is initialized with a realistic salinity-temperature-depth (STD) 
profile acquired in the northern Gulf of Mexico at the Mobile OTEC site in the 
summer of 1977 (Starr and Maul, 1982). The maximum Brunt-Vaisala frequency is 
roughly 10 cph (Fig. 6.1) at 50 m, then decreases to roughly 2 cph over a vertical scale 
of 250 m and is a constant 1.5 cph below 400 m. This initial profile closely resembles 
the Brunt-Vaisala frequency profile estimated from climatology and the AXBT 
observations of Black (1983) (see Fig. 2.4b), except that the maximum Brunt-Vaisala 
frequency is shifted downward by 10 m. 

The first five free baroclinic modes are calculated using the Brunt-Vaisala 
frequency profile and the Sturm- Li ouville problem (Fig. 6.2). Rapid changes in the 
vertical structure of the <p n ^s occur below the mixed layer for modes 3, 4 and 5. The 
corresponding phase speeds and horizontal wave lengths are quite similar to those 
calculated for the Brunt-Vaisala frequency profiles derived from the AXBT's (see Fig. 
2.4). For example, the phase speeds are 2.9, 1.3, 0.8, 0.5 and 0.4 m/s for the first five 
baroclinic modes for the OTEC profile. 

B. SIMULATED CURRENT AND TEMPERATURE STRUCTURE 

The pattern of simulated mixed layer velocities (Fig. 6.3) is similar to that 
predicted from linear theory (Geisler, 1970). Upwelling and downwelling are associated 
with the convergent and divergent velocities. The wavelength of these cells is given by 



entrainment heat flux in a u* 3 formulation. 
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Figure 6.1 Brunt-Vaisala frequency (cph) profile from an OTEC STD station 10 at 29 
°N and 88 °W. 

product of the IP and U h . The effect of the nonlinearities is to displace the velocity 
maximum to the right of the storm track (Chang and Anthes, 1978), since linear 
simulations using the same Rankine vortex have the maximum velocities along the 
storm track. Maximum velocities of 140 - 160 cm/s are simulated in the directly forced 
region. At the same distance (r = 3-4 R ) from the storm center as Frederic data 
were acquired, simulated current velocities range between 100- 120 cm/s, which agrees 
with the observed mixed layer currents at CM A3 (see Fig. 2.4). 

The pattern of velocities in the thermociine (285 m) differs considerably from the 
velocities in the upper ocean (Fig. 6.4). Maximum velocities are about 16 cm/s and are 
towards the west. However, the eastward velocities (located on the left side of the 
track) are elongated in the cross-track direction. These velocity perturbations are in 
(out) of phase with the baroclinic ridge to the right (left) of the storm track. 

1. Along-Track Sections 

Along-track sections of the v-component of velocity and temperature at 
x = 4R are shown in Fig. 6.5a. The vertical coordinate (l,) is scaled using WKBJ 
theory (3.7) to emphasize the processes occurring in the upper 200 m. This scaling 
stretches the grid spacing even more in the thermociine and shrinks it in the bottom 
layers. The horizontal wavelength in the simulations corresponds to the linear theory 
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Figure 6.2 Amplitude of the horizontal velocity eigenfunctions for : mode 1 (solid), 
mode 2 (dashed), mode 3 (dotted), mode 4 (chain-dotted) and mode 5 (chain-dashed) 
based on Brunt-Vaisala frequency (Fig. 6.1). 

wavelength of 580 km. Within the first IP (x= 2A), the ocean currents are excited 
throughout the water column because the free surface slope induces a depth- 
independent velocity component. This is the component that previous numerical 
modeling studies have ignored by imposing a rigid lid in the model. Notice that the 
slope of the phase lines is upward, which represents upward phase propagation and a 
downward energy propagation (Leaman and Sanford, 1976). Maximum velocites in the 
thermocline increase to 20-30 cm/s near 5 IP. The u-component also shows similar 
behavior in the phase tilt and amplitudes of the currents. 

The along-section temperature changes (Fig. 6.5b) are modulated by the near- 
inertial cycle in the thermocline with maximum changes are about 1-1.2 °C. There is a 
rapid decrease in the mixed layer temperature during storm passage. The temperature 
increases in the thermocline are associated with down welling processes outside of the 
main upwelling zone along the storm track. 
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Figure 6.3 Simulated mixed layer u-component of current (cm/s) from the nonlinear 
model with a contour interval of 20 cm's. The axes are scaled in terms of R max (30 
km) and A (580 km) and the storm is moving westward along y = 0 at 6.5 nvs. The 
storm center is located at x = A. 

2. Effects of the Free Surface 

The deflection of the free surface simulated by the model is shown in Fig. 6.6. 
The maximum changes in the free-surface height are 12-14 cm and are centered along 
the storm track. This deformation in the free surface is caused by the mean mass 
divergence in the upper ocean which subsequently changes the mean vertical pressure 
gradient force. Thus, currents are excited throughout the water column via the 
pressure gradient force. 

To understand the contribution of the free surface slope to the simulated 
ocean currents, the geostrophic velocity components are calculated from the gradient 
of the free surface (Fig. 6.7). On the right (left) side of the track, the maximum u- 
components arc 12-14 envs in the direction of (opposite to) storm movement. These 
velocities are similar to those observed in the thermocline (Fig. 6.4). There also 
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Figure 6.4 Simulated thermocline (285 m) u-component of current (cm/s) from the 
nonlinear model with a contour interval of 4 cm/s. The axes are scaled as in Fig. 6.3. 

appears to be small cells within y= =k4R max that have a wavelength of 580 km. The 
maximum v-component velocities (Fig. 6.7b) are slightly weaker with amplitudes of 4-8 
cm/s. However, there are also oscillations in the v-component that have near-inertial 
wavelengths of 580 km. Thus, near-inertial oscillations are included in the barotropic 
response. 

In summary, these simulations have the rightward bias of the maximum 
current response at about r = 2-4 with maximum mixed layer currents of • 

120-140 cm/s in that region. Furthermore, the near-inertial wavelengths are about 580 
km as predicted by linear theory. Thus, the simulated currents appear to have 
sufficient veracity to allow a comparison to the moored current meter observations. 
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Figure 6.5 (a) Simulated v-component of current (cm/s) and (b) AT (°C) from the 
nonlinear model at x = 4R max . The abscissa is scaled as in Fig. 6.3. and the ordinate is 
WKBJ-stretchcd using (3.7). 
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Figure 6.6 Deviations of the free surface height (cm) simulated by the model relative to 
the undisturbed height. The axes are scaled as in Fig. 6.3. 

C. DATA ANALYSIS 

In this section, the simulated velocities are analyzed to determine the frequency 
of the near-inertial response. The frequencies of the waves should be blue-shifted 
above the local inertial frequency by the amount predicted from the mixed layer and 
thermocline Burger numbers (Price, 1984). In the Frederic observations, the vorticity 
of the mean flow could not be calculated because only two or three moorings were 
available. However, the mean flow vorticity can be estimated from the numerical 
simulations and compared to the diagnosed frequencies of the forced waves. 

To simplify comparisons of the numerical simulations with the observed current 
meter observations in Frederic, the simulated data are rotated into a coordinate system 
with the storm moving towards the north. This can be done by simply rotating the 
simulated velocities clockwise by about 90 °. Since the model assumes an f-plane, 
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Figure 6.7 (a) The u-component (cm's) and (b) the v-component (cm/s) of current 
induced by the slope of the free-surface simulated by the model. The axes arc scaled as 
in Fig. 6.3. 
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beta-dispersion of inertial waves is not possible in the numerical simulation (Anderson 
and Gill, 1979). 

1. Frequency Analysis 

The horizontal velocities are first least-squares fitted with (4.1) to the model of 
Mayer et al. (1981) to diagnose the frequency of the oscillations over a 5 IP segment. 
The variance reduction curves from the simulated model velocities along x = 4R max are 
shown in Fig. 6.8. These variance reduction curves are very similar to those in Fig. 4.4, 
especially in the mixed layer. As noted in Chapter IV, the carrier frequency may not 
necessarily be aligned with these maxima, since it is defined to be the frequency that 
minimizes the covariance of the residual signals. Based on the mixed layer Burger 
number using the model parameters, the predicted frequency shift at 25 m is about 0.02 
to 0.03f. The least-squares fit model can account for about 80% of the variance 
between 1.02 to 1.03 f and over 60 % of the simulated current variability in the upper 
part of the thermocline at 105 m. At 175 m (not shown), the maximum variances 
accounted for by the model are only about 45 %. Note in Fig. 6.8b that the horizontal 
velocity eigenfunction for the first mode is nearly zero between 135 and 175 m. At 
deeper levels (Fig. 6.8c and d), the least-squares model can again account for more of 
the variance of the simulated currents. 

The diagnosed frequencies are generally blue-shifted between 1-6% above the 
local inertial frequency (Table 16), which is consistent with the mixed layer Burger 
number. This blue-shift in the frequency is evident in the upper ocean as the carrier 
frequency slowly changes from 1.01 fat the surface (5 m) to about 1.04 - 1.05 f at 135 
m. Below the mixed layer, Price (1984) suggests that there is a blue-shift with depth 
that is proportional to the thermocline Burger number. Kundu and Thomson (1985) 
also found a consistent blue-shift with depth in the upper 200 m. The thermocline 
Burger number in the model is about 0.36 and the corresponding frequency shifts are 
proportional to T/4 or roughly 0.09f. Even with the uncertainities in the least-squares 
model of 0.02f, the frequency shift at 135 m is only 1.06f. Moreover, below this 
maximum, the diagnosed frequencies decrease slightly with depth, rather than increase 
as predicted from the scaling arguments. If only 3 or 4 selected levels are available as 
in Price (1984), a consistent blue-shift with depth might be inferred. However, it is 
clear the diagnosed frequencies do not show a consistent blue-shift with depth either in 
the simulated currents (see Table 17) or in the Frederic observations. 
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Figure 6.8 Explained variance (%) for the simulated cross-track (solid) along-track 
(dashed) velocity components at x = 4R max For trial Frequencies normalized by the 
local inertial Frequency (<r/F) at a) 25 m, b) 105 m, c) 285 m and d) 435 m. 
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TABLE 16 



Near-inertial frequencies (periods) at which a maximum 
response occurs in the simulated currents based on a least- squares fit over 5 



Depth 

(mj 


IP. 

*/f 


Period 

(h) 


5 


1.01 


24.5 


15 


1.02 


24.2 


25 


1.02 


24.2 


35 


1.02 


24.2 


45 


1.02 


24.2 


55 


1.03 


23.9 


70 


1.03 


23.9 


90 


1.03 


23.9 


105 


1.04 


23.7 


135 


1.04 


23.7 


175 


1.02 

1.03 


24.2 


225 


23.9 


285 


1.03 


23.9 


355 


1.03 


23.9 


435 


1.04 


23.7 


515 


1.06 


23.3 


595 


1.06 


23.3 



2. Mean Flow 

One possible mechanism for shifting the frequency of the near-inertial waves is 
the superposition of a sheared mean flow (Mooers, 1975; Kunze, 1985). If the 
variations of the mean flow are comparable to the wavelength of the propagating near- 
inertial waves, the relative vorticity of the mean flow can effectively shift the frequency 
of the near-inertial waves above (below) f in regions of cyclonic (anticyclonic) vorticity. 

The simulated current fields are averaged using a simple running mean over an 
inertial wavelength (Fig. 6.9a). The cross-track scale of the mean flow in the mixed 
layer induced by the hurricane is about 8-10 R max or about 300 km. The maximum v- 
component of velocity is 30-40 cm/s under the storm forcing. On the left side of the 
track, the v-component flows southward at about 10 cm/s whereas the flow is northerly 
at 10 cm/s on the right side of the track. In the mean currents at CM A3 in Frederic 
(see Fig. 4.3a), the v-component initially changed by 8-10 cm/s and then reversed 
direction towards the south within an IP to a maximum of 10-15 cm/s. In the model 
simulations, the magnitudes after an IP (y = -2A) are comparable, but are flowing in 
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opposite directions. By contrast, the mean u-component has a similar cross-track 
scale, but damps within an IP. Although some time-mean flows are evident in the deep 
layers, the magnitudes of the flow are only about 2-4 cm/s and the cross-track scales 
are about 2-4 R max - 

The relative vorticity (Q of the mean flow is calculated and normalizied by the 
local Coriolis parameter (2f) to form an effective Coriolis frequency shift (Kunze, 1985) 




In the mixed layer (Fig. 6.9b), the cross-track width of the f e fp is only 120 km, as 
compared to the scale of the near-inertial oscillations of 580 km. Kunze (1985) shows 
that the scale of the vorticity field has to be comparable to the wavelength of the 
oscillations for the vorticity to shift the observed near-inertial frequency. The f £ |p has 
a maximum value of about 0.02-0.03 f along the track and is cyclonically rotating 
w’hich can also blue-shift the frequency of the waves. There are two small cells of 
anticvclonically rotating vorticity on either side of the track. However, they are 
transient features and damp quickly. As noted in Chapter IV, the vorticity of the 
mean flow could not be resolved with only the current meter arrays in Frederic. 

The estimated frequency shifts from Table 17 and the f e |y are compared over 
depth in Fig. 6.10. It is clear that at x = 0, 2R max and 4R max , the shifts in cr/f near 
the surface are out of phase with the f e ff vertical structure. In the upper levels (5-15 
m), the f e {pis a maximum of 0.07 f, whereas the diagnosed frequencies are only 
1.01-1.02 f. This estimate of f e fymay be somewhat misleading since the upper levels 
(5-15 m) are in the diffusion dominated zone in the model. As depth increases, f e pp 
decreases to nearly zero below 55 m. By contrast, the cr/f increases with depth to a 
maximum of 1.08 at 135-175 m along the track (x= 0). A similar increase with depth is 
apparent at x = 2R max . However, the maximum cr/f at 4R max is only 1.04. The 
nearly uniform increase with depth of the cr/f agrees with predictions calculated from 
the Burger numbers (Price, 1984). The frequencies of the simulated near-inertial 
oscillations decrease with depth even directly along the track. Thus, the consistent 
blue-shift with depth may be limited to the upper ocean with vertical scales comparable 
to the thermocline. 

D. NEAR-INERTIAL RESPONSE 

The simulated velocity data are demodulated at the carrier frequencies given in 
Table 17 to isolate the near- inertial response in the simulated current fields. The 
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nondimcnsionally sealed in terms of R max (30 km) and A (580 km) with the storm 
center now at (0,-1). 




o/f C & 

0 05 1 1 05 1.10 - 0.05 0.00 0.05 0.10 




Figure 6.10 Vertical profiles of the normalized diagnosed frequencies <J, f (left) and the 
f e ^ (right) at x = 0 (solid), x=2R max (dotted) and x = 4R max (dashed) based on 
simulated velocities from the 17-level model. 

technique is similar to that used in Section 4.2 for the observations, except that the 
data are filtered using a simple running inertial period, rather than the Lanczos filter 
used in the observations. This change in the technique is necessary because the 
simulated time series only extends for 5 IP. 

1. CW-rotating amplitudes 

Since the demodulation process yields four amplitudes (cosine and sine for 
each velocity component), the CW-rotating amplitudes (Fig. 6.11) can be obtained 
following the treatment by Gonella (1972) and Mooers (1973). As expected, the largest 
CW-rotating amplitudes are to the right of the track in the upper ocean and exceed 100 
cm/s in the mixed layer (Fig. 6.1 la). The amplitudes decrease over e-folding scales of 
2-3 IP, which agrees with the observations (see Fig. 4.5). For example, the maximum 
amplitudes in the mixed layer at CMA3 were about 85 cm/s after removal of the 
surface wave velocity components. Thus, there is fairly good agreement between the 
CW-rotating amplitudes in the upper ocean between the observations and the model. 
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At 70 m (Fig. 6.11b), the maximum CW-rotating amplitudes are about 50-60 cm/s, or 
about half the amplitude in the mixed layer. There is also a shift of about an IP 
between the maximum CW-rotating amplitudes in the mixed layer and thermocline. 

This time shift suggests a phase tilt with depth that is expected for vertical energy’ 
propagagtion. These amplitudes substantially agree with the amplitudes of 18-22 cm/s 
observed at 250 and 440 meters at CMA3. In the deep layers (Fig. 6.1 lc and d), there 
is a fairly signficant amplitude of 16-18 cm/s in the near-inertial frequency band! 
Presumably, these large amplitudes are due to the free-surface boundary condition. 

The CCW-rotating amplitudes are 2-8 cra/s (not shown) and when combined with the 
CCW amplitudes form near-inertial ellipses that enable energy to propagate vertically 
(LeBlond and Mysak, 1978). 

The four amplitudes are combined to form the characteristics of the CW- 
rotating amplitudes and phases for direct comparison to the observations. In the 
mixed layer, the amplitudes and phases are fairly close. However, the simulated 
amplitude in the thermocline (180 m) is only 4 cm/s as compared to 15 cm/s observed 
at CMA2. There is a node point in the mode 1 amplitude (see Fig. 6.1) w’hich may 
account for a large percentage of the difference. The observations and model near- 
inertial currents do seem to be in phase at 1 80 m. At 285 m, the model amplitudes are 
about 5 cm/s less than those observed and there is a difference in phase of about 100 °. 
How r ever, the amplitudes at 435 m are underestimated, but the phases are within 2 °. 
Some of the discrepancies at these two levels have to be attributed to the sloping 
bottom, although the phases at the bottom are nearly equal. Although the model 
simulations have slightly less amplitude than the observations, it is encouraging that 
the model does predict significant currents below’ 200 m! 

2. Modal Analyses 

The simulated velocities in the near-inertial wave band are then least-squares 
fit (Fig. 6.12) to the barotropic mode (not shown) and first three baroclinic modes in 
Fig. 6.2. Notice that there is indeed a barotropic component of 4-6 cm/s within the 
near-inertial band, which agrees fairly w’ell w’ith the depth-averaged component 
estimated from the moored current meter arrays of 7-10 cm/s (see Fig. 4.12). 

Moreover, the mode 1 amplitudes along x = 4R max (Fig. 6.12b) is a maximum of about 
30 cm/s as compared to 24 cm/s in the observations. The second baroclinic mode 
amplitudes (Fig. 6.12c) are also within 3-4 cm/s of the those estimated from the fits to 
the observational analyses. Unfortunately, the mode 3 amplitude could not be resolved 



127 



-y/A or t/IF -y/A or t/!P 





F igure 6.11 CW-rotating amplitudes from the simulated, demodulated velocities (cm/s) 
for the a) 25 m, b) 70 m, c) 285 m and d) 435 m. The axes are scaled as in Fig. 6.9 
with storm center at (0,-1). 
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TABLE 17 



Comparisons of the CW-rotating amplitudes and phases 
between the Frederic observations and model simulations based on a 5 IP 

average. 



Observed Model 



Depth 


CW 


O 


CW 


O 


(m) 


(cm/s) 


(°) 


(cm/s) 


(°) 


25 


70 


' 142 


60 


104 


175 


15 


-103' 


4 


-131 


285 


20 


27 


14 


-74 


435 


19 


-69 


14 


-71 



with the limited vertical sampling by the current meter moorings. However, the model 
suggests that the third baroclinic mode amplitudes are 4-6 cm/s. 

E. SUMMARY 

The simulated near-inertial response to a hurricane-like wind stress has very 
similar characteristics to the observations acquired in Frederic. For example, there is a 
blue-shift in the frequencies of the waves of 1-8% above f. It appears that the vorticity 
of the mean flow is not responsible for the frequency shifts. The mean (time) flows set 
up by the hurricane have similar magnitudes in the model and observations. The CW- 
rotating amplitudes throughout the water column also have similar magnitudes which 
is especially apparent in the near-bottom amplitudes of 16-20 cm/s. 

The free surface slope that is induced by mean divergence of the currents excites 
a depth-independent component in the near-inertial frequency band. This barotropic 
component has amplitudes of 4-6 cm/s compared to 7-10 cm/s derived from the 
observations. Furthermore, the amplitudes of the first two baroclinic modes are within 
3-5 cm/s of those estimated from the least-squares fits. Thus, it appears that the free 
surface model is simulating realistic ocean currents in response to strong atmospheric 
forcing. 
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Figure 6.12 Velocity coefficients from the simulated, demodulated velocities (cm/s) for 
the a) barotropic mode, b) mode 1, c) mode 2 and d) mode 3. The axes are scaled as in 
Fig. 6.9 with storm center at (0,-1). 
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VII. CONCLUSIONS 



The ocean current data collected during and subsequent to the passage of 
hurricane Frederic represented a unique opportunity to study the time evolution of the 
forced three-dimensional response. The Frederic data set included sufficient spatial and 
temporal sampling to understand the low-mode (large vertical wavelength) near-inertial 
response. A striking feature was the fairly good agreement in the near-inertial 
amplitudes with the time scales for the first two baroclinic modes as predicted from a 
linear, inviscid model (Gill, 1984). For instance, the first baroclinic mode maxima 
occurred within an IP and decreases rapidly as energy is propagated out of the source 
region in the mixed layer. A secondary increase in mode 1 amplitude at about 3-4 IP 
corresponds to the predicted time for the first mode to become in-phase with the higher 
order modes. This secondary increase in the mixed layer amplitudes about 3-4 IP 
following storm passage at the arrays also agrees with the numerical simulations of 
Kundu and Thomson (1985). Secondary increases in the clockwise rotating, near- 
inertial amplitudes 7-8 IP after storm passage were associated with the second 
baroclinic mode in the thermocline and near-bottom layers. Amplitudes in the 
thermocline started to decrease substantially after 15 IP which is consistent with 
predictions from Gill's theory. 

The explained variances associated with the first and second baroclinic modes 
were nearly equal at Current Meter Array CMA3 when averaged over 7-10 IP, whereas 
the second mode dominated the response at CMA2 by almost a factor of three. These 
differences can be explained in terms of the differences in the bottom topography at the 
arrays. The inclusion of the bottom slope in the normal mode decomposition of the 
current meter data affects of the baroclinic modes that can be resolved in the current 
meter data. The first baroclinic slope-mode explained about 45 % of the near-inertial 
variance or roughly 5 % more than the first baroclinic, flat -bottom mode. The amount 
of variance contributed by the second baroclinic mode was about 6 % above the 
variance explained in the flat-bottom case. These increases in the explained variances 
due to the inclusion of slope effects generally agree with those by Lai and Sanford 
(1986). 



131 



One of the significant features observed in all of the current meter measurements 
was a rapid current increase throughout the column. Bottom currents increased 
significantly in the DeSoto Canyon region and at the OTEC site within hours of the 
upper ocean response, which began roughly 9 h prior to the point of closest approach 
of hurricane Frederic in the DeSoto Canyon region. Brooks (1983) also observed an 
ocean response to 700 m during the passage of hurricane Allen. Thus, the Frederic 
data provide important confirmation that a deep ocean response occurs much more 
rapidly than had been previously predicted on theoretical grounds. 

Perhaps the most controversial aspect of the Frederic moored current meter data 
is the appearance of the depth-averaged component. Because of limited vertical 
sampling, these estimates of the depth-averaged response to hurricane Frederic may be 
contaminated by aliasing of higher-order baroclinic modes. Nevertheless, the depth- 
averaged flow accounted for about 20 % of the near-inertial variance at all the sites 
over the first 10 IP. The initial increase in the depth-averaged current (8 cm/s) may 
have been associated with a free-surface slope since the current increases were detected 
to the bottom at all arrays. However, these initial increases in the depth-averaged 
current only account for 30 % of the observed variance. It is concluded that the 
depth-averaged currents of 7-10 cm/s in the Frederic moored current meter arrays are 
consistently too large to be an aliasing problem, and thus are likely to be a real 
physical aspect of the ocean response to hurricanes. 

The linear, two-layer model of Geisler (1970) has been extended to include a 
continuously stratified fluid with the wind stress superposed onto the first five 
baroclinic modes (Kundu and Thomson, 1985). Expressions for the 3-dimensional 
velocity structure have been derived that include the wind stress curl, the stress 
divergence and Ekman-type terms. The model hurricane is a Rankine vortex based on 
observed storm parameters. The forcing functions for each mode are convolved with 
the Bessel function to determine the horizontal structure of the 3-dimensional velocity 
field. In general, the wind stress curl dominates the ocean current response. However, 
the Ekman and stress divergence terms that were not treated in Geisler (1970) 
contribute to the near-inertial response within the direct forcing regime. These terms 
dampen quickly within one or two wavelengths of the forced regime, whereas the curl 
has more of a persistent effect on the ocean response. 

For the first time, the vertical structure of the ocean current and temperature 
response to the passage of major hurricanes has been mapped using AXCP 
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measurements (Sanford et ai, 1987). These experiments demonstrate that AXCP's can 
provide ocean current and temperature data under extreme conditions. In this 
research, the vertical structure of the near-inertial response to hurricane Norbert is 
isolated from the ocean current profiles by first subtracting a spatially-averaged mean. 
The near-inertial profiles are then fit to the first three or four forced baroclinic modes 
from the extension of Geisler's model. The forced, linear model explains over 70 % of 
the observed current variance in the right-rear quadrant of the storm with only the first 
four baroclinic modes. Near-inertial processes appear to contribute to an enhanced 
vertical shear feature observed near the base of the mixed layer in most of the profiles 
under the direct forcing of the hurricane. In terms of the energetics, the KE based on 
the simulated profiles represent the observed KE reasonably well (67-80 %) on the left 
side of the storm. Thus, the major portion of the vertical structure of the ocean 
current response can be represented by only three or four vertical modes. Mixing and 
nonlinear effects that have been neglected in this treatment may be important in some 
areas. 

Hopkins (1982) had attempted to simulate the ocean response to hurricane 
Frederic with a rigid lid, baroclinic model. Although there was good agreement with 
the mixed layer currents, the currents in the thermocline were underestimated in 
comparison to the observations. To assess the relative role of the barotropic (depth- 
averaged) and baroclinic modes on the oceanic response to hurricane, the near-inertial 
response to a hurricane- like wind stress is isolated from numerical simulations from a 
free surface model (Chang, 1985). This is the first model that includes both the 
barotropic and baroclinic response to hurricane forcing. The maximum simulated 
currents are displaced between r = 2-4R max to the right of the storm track, which 
agrees with previous numerical studies. The wavelength of the forced waves is equal to 
the predictions from linear theory. Thus, the model simulations appear to have 
sufficient veracity to allow comparisons to observations. In the free surface model 
simulations, the near-inertial response is detected throughout the water column because 
the slope of the free surface induces a depth-independent response. Typical bottom 
amplitudes of the CW-rotating currents were 15 cm/s compared to 18 cm/s observed in 
Frederic. Using the same techniques used with the observations, the near- inertial 
amplitudes are fit to the first few vertical modes. The modal amplitude of the 
barotropic component is 4-8 cm/s as compared to the estimate of 7-10 cm/s derived 
from the observations. Thus, it is suggested that the simulated barotropic response 
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supports the observations of the depth-averaged response despite the limited vertical 
sampling in the moored current meter arrays. There is also considerable agreement 
with the maximum amplitudes of the first tw r o baroclinic modes in the simulations. 

The combination of moored current meters to provide the detailed time evolution 
and the AXCP's to describe the instantaneous vertical structure in great detail has been 
very effective in increasing our understanding of the ocean response to hurricane 
forcing. It is particularly satisfying that the theoretical treatments agree in so many 
aspects with the observed time evolution and the vertical structure, including (for the 
first time) the forced response under hurricane Norbert. Simulations with realistic 
forcing provide a more complete three-dimensional and time-dependent representations 
of the ocean response. Applying the same analysis techniques as used in the 
observational data results in a satisfactory confirmation of many aspects of the ocean 
response. Thus, we have considerable confidence that a correct representation of the 
vertical structure of the ocean response has been derived from the sparse moored 
current meter and AXCP observations. 
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APPENDIX A 

THREE-DIMENSIONAL VELOCITIES 



1. VERTICAL VELOCITY 

The vertical velocity horizontal structure function will be treated first because of 

it's implications with upwelling and downwelling regimes in the wake of a moving 

✓ 

tropical cyclone. Differentiating (3.10) with respect to time and multiplying (3.11) by 
the Coriolis parameter f and adding yields: 
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The corresponding equation for the v-component is given by: 
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After differentiating and combining the terms, and dividing by c n 2 the expression 
becomes 
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where V 2 is the horizontal Laplacian given by 
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and Vxx n = dYJdx - dXJdy and V*T n = dXJdx +dYJdy which represent the wind 
stress curl and divergence, respectively. Substituting for the pressure term from (3.15), 
the governing equation for the vertical velocity for the n** 1 baroclinic mode becomes 



a 2 a 2 

dx 2 dy 2 



1 

c 



a 2 



n 



2 at 2 



2 >W = fVxl H V *T 



n 



at 



n 



(A.l) 



135 



where Vxt n = dY n ;dx - dXJdy and V*r n = dXJdx + dYJdy . Equation (A.l) is 
analogous to equation (10) in Geisler (1970) except that the phase speed (c n ) and the 
inverse deformation radius (a ) are for the n 1 * 1 mode. 



2. HORIZONTAL VELOCITY 

The spatial distribution of the horizontal velocity eigenfunctions is sought using 
(3. 10), (3. 11) and (3.15) (Longuet-Higgins, 1965). Differentiate equations (3.10 and 
3.11) with respect to time and substitute for pressure from (3.15): 
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A vorticity equation is formed by cross differentiating equations (A. 2) and (A.3) 
dt 2 

and solving for v n yields 
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Next, take the partial with respect to time of (A. 2) and add it to the product of 
Coriolis parameter and (A.3) to form 
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and taking the partial once more with respect to time and substituting (A.4) for v 
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Dividing by c ^ and integrating twice with respect to time forms 
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Equation A. 7 can then be simplified further invoking the steady-state and transforming 
into the coordiate system defined in Chapter III. Thus, the governing equation for the 
horizontal structure of the u-componcnt is 
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(A. 8) 



and the governing equation for the v-component is developed in a similar manner 



In these equations, the second terms represent those components that are driven by the 
wind stress, Y /f and -X /f or the Ekman type velocities. The third terms in the u and 
v equations represent a sink of vorticity and the y-component of the divergence, 
respectively. Similarly, the effect of the divergence term is to augment the curl term 
because the minus sign makes the term of the same sign. 




.(A.9) 
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APPENDIX B 

HYPERBOLIC OPERATOR 



The inversion of the integral (page 53) yields a fundamental solution or a Green's 
function (Stakold, 1979). The double integral (INT) is given by 



I NT = 
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as defmed in the Chapter III above. This integral can be recast into the cosine 
transforms 
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The transform in brackets is given by 
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(Erdelyi et al 1954) which renders the integral equation as 

IXIT 1 « sin(V(k 2 +l)y') „ 

INT = — | — = cos(kx ) dk . 

271 l (k 2 + 1) 

This inverse transform can be found in Erdelyi et al., (1954) as 



“ J 0 ( (y' 2 -x ,2 ) 1/2 ) 

where J Q represents the zero order Bessel function. Thus the inversion of the transform 
yields 

INT= ^ . 



The Bessel function of order 0 is the fundamental solution (Green's function) to the 
integral equation which is valid for 0 < x' < y\ In terms of the atmospheric scaling 
(2R max ), the inverse transform is 

INT = ((y' 2 -x' 2 ) 1/2 ) . 
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APPENDIX C 
ERROR ANALYSIS 



The amount of aliasing of low-mode, baroclinic energy by limited vertical 
sampling is addressed using the predicted currents from a 15-level, primitive equation 
model (Adamec et ai, 1981). The model is externally forced by a fast-moving 
hurricane (8.5 m/s) and integrated over 3 IP. Since the storm speed exceeded the first 
baroclinic mode phase speed (3 m/s), a baroclinic response is expected (Geisler, 1970). 

A rigid lid is imposed in the model and the vertically-averaged flow is required to be 
zero at each grid point. The variance reduction method was performed on the depth- 
time series over 3 IP in the near-inertial frequency band. A depth series containing the 
real part of the velocity was formed prior to the modal decomposition. 

To assess the amount of bias due to limited vertical sampling, the depth series 
was vertically averaged using all the 15 and only 3 levels of data respectively. Using 15 
levels, the depth-averaged, near-inertial component was about 3 cm/s. Thus, this 
depth-averaged component is assumed to be correct since the model requirement of a 
zero vertical average is valid only over the total frequency range. The depth series was 
subsampled at levels that correspond to the current meters on the NAVOCEANO 
arrays (3 levels). The depth-averaged, near-inertial component was 1 cm/s or a change 
of 2 cm/s from the 15 level case. Thus, limited vertical sampling induces a 2 cm/s bias 
in the estimation of the depth-averaged current. 

A depth-averaged component is subtracted from the near-inertial profile prior to 
least-squares fitting of the first two baroclinic modes. The first mode coefficient in the 
15-level case was 17 cm/s as compared to 23 cm/s for the 3-level experiment. Thus, the 
bias in the first mode amplitude due to limited vertical sampling is about 6 cm/s. The 
second baroclinic mode amplitudes were 14 and 10 cm/s for the 15 and 3-level 
experiments, respectively. Thus, the bias in amplitude with only 3 observation levels is 
approximately 4 cm/s for the second baroclinic mode. It should also be noted that first 
baroclinic mode was overestimated and the second mode coefficient was 
underestimated due to limited vertical sampling. The total amount of energy in the 
vertical average and first two modes was conserved in both experiments. In summary, 
the biases of this magnitude in the depth-averaged and first two baroclinic mode 
amplitudes are not large enough to negate the conclusions of this study. It is 
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concluded that Frederic excited a non-zero depth-averaged component and energetic 
baroclinic modes 1 and 2 that have been inferred from the limited ocean current data. 
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